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Abstract 

Interest in short waves on the ocean surface has been growing over the last three decades because 
they play an important role in surface electromagnetic (e.m.) scattering. Currently radars and scat- 
terometers which use e.m. scattering to remotely examine the ocean can produce estimates of the 
surface wind field, surface currents, and other scientifically important ocean processes. These esti¬ 
mates are based on models which depend on a thorough understanding of electromagnetic scattering 
mechanisms, and of the three-dimensional surface wave field. Electromagnetic scattering theory is 
well developed, but the short wavelength portion of the surface wave field has only recently been 
experimentally explored. A single, consistent, and accurate model of the energy distribution on the 
ocean surface, also known as the wave height spectrum, has yet to be developed. 

A new instrument was developed to measure the height of waves with 2-30 cm wavelengths at an 
array of locations which can be post-processed to generate an estimate of the two-dimensional wave 
height spectrum. This instrument (a circular wire wave gage buoy) was deployed in an experiment 
which gathered not only in situ measurements of the two-dimensional wave height spectrum, but 
also coincident scatterometer measurements, allowing the comparison of current e.m. scattering and 
surface wave height spectrum models with at sea data. 

The experiment was conducted at the Buzzards Bay Tower located at the mouth of Buzzards 
Bay in Massachusetts. A rotating X-band scatterometer, a sonic anemometer, and a capacitive wire 
wave gage were mounted on the tower. The wave gage buoy was deployed nearby. The resulting data 
supports a narrowing trend in the two-dimensional spectral width as a function of wavenumber. Two 
current spectral models support this to some extent, while other models do not. The data also shows 
a similar azimuthal width for the scatterometer return and the width of the short wavelength portion 
of the wave height spectrum after it has been averaged and extrapolated out to the appropriate Bragg 
wavelength. This appears to support current e.m. composite surface (two-scale) theories which 
suggest that the scattered return from the ocean at intermediate incidence angles is dominated by 
Bragg scattering which depends principally on the magnitude and shape of the two-dimensional 
wave height spectrum. However, the mean wind direction (which corresponds well with the peak of 
the scatterometer energy distribution) and the peak of 20 minute averages of the azimuthal energy 
distribution were out of alignment in two out of three data sets, once was by nearly 90°. There are a 
number of tenable explanations for this including instrument physical limitations and the possibility 
of significant surface currents, but none that would explain such a significant variation. Given that 
there are so few measurements of short wave directional spectra, however, these results should be 
considered preliminary in the field and more extensive measurements are required to fully understand 
the angular distribution of short wave energy and the parameters upon which it depends. 
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Chapter 1 


Introduction 


The emergence of airborne and space borne remote sensing capabilities have provided a long- 
awaited tool for studying the surface and sub-surface properties of the world’s oceans. This 
is because electromagnetic backscatter from the sea surface depends chiefly on the ocean 
surface roughness at both large and very small scales. The surface roughness is affected 
predominantly by the wind field, but also by other factors including surface currents, ocean 
floor topography, and surface films. Electromagnetic scattering techniques operated on 
satellite or airborne platforms have the capability, therefore, of providing information about 
these ocean processes on potentially synoptic scales. Understanding the mechanisms by 
which electromagnetic measurements are coupled to the surface roughness via the primary 
mechanism of Bragg scattering, and through the surface roughness to the actual wind field, 
surface currents, etc., however, is a very difficult task which requires detailed information 
about the surface wave field. While the large gravity waves which contribute to surface 
roughness have been under investigation for many years, technological advances allowed 
the study of the smaller capillary and capillary-gravity waves to begin in earnest only in 
the last three or four decades. As the theories describing electromagnetic wave interactions 
with the ocean surface have improved, however, so has the need for more and more detailed 
spatiotemporal information about the surface wave field. 

The directional spreading of the surface spectrum has historically been the weakest 
aspect of spectral and scattering models because of the difficulty in obtaining at-sea mea¬ 
surements of it at small enough wavelengths. Toward that end, a new instrument has been 
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developed to aid in the study of the directional characteristics of these very sort waves on 
the ocean surface. It is known as the circular wire wave gage. It is an autonomous buoy 
capable of providing directional energy distribution information about waves with wave¬ 
lengths between 2 and 30 centimeters. This is the wavelength range which impacts the 
Bragg scattering mechanisms important in microwave backscatter. It was deployed on Oc¬ 
tober 28, 1998 at the mouth of the Buzzards Bay in Massachusetts as part of an experiment 
designed to compare in situ measurements of both the surface wave height spectrum and 
electromagnetic backscatter with currently proposed models for both of these processes. 

In this work the current theoretical knowledge about the surface wave height spectrum 
and electromagnetic backscatter from the ocean surface is coupled with in situ experimental 
data, providing evidence not previously available about the applicability and limitations of 
these theories. Chapter 2 presents the current surface wave and electromagnetic theories 
as well as describing other available techniques and instruments which have been used to 
measure short surface waves. A detailed description of the new instrument being used to 
measure directional distribution of short waves is given in Chapter 3. Chapter 4 describes 
the instruments deployed and the data collected during the Buzzards Bay Tower experiment. 
And finally, Chapter 5 ties the theory to the experimental data to illustrate the strengths 
and weaknesses of both. 
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Chapter 2 


Background and Theory 


Accurate analytical and predictive models for ocean surface waves have been under develop¬ 
ment for many years. Long waves have historically received most of the attention because of 
the large amount of energy they contain, as is illustrated in Figure 2-1. The shorter waves 
saw little attention until World War II when the Navy became interested in understanding 
’sea clutter’, a radar phenomenon which could mask radar signatures of surface ships and 
submarines caused by wind roughened water surface [14, 69]. Since that time much has 
been learned about short waves and their influence on important fields of study such as 
air-sea interactions (e.g. gas transfer) and wave growth and generation. 

Many in situ and remote instruments have been developed to aid the study of these 
short waves. In order to understand the strengths and limitations of these instruments it 
is necessary to appreciate the physical phenomena they are measuring. The first section 
of this chapter, therefore, describes surface wave propagation and energy distribution. The 
remote sensing techniques used to study ocean surface waves require knowledge of the 
interaction of electromagnetic waves with the water surface. Accordingly, the second section 
describes electromagnetic scattering and radar cross section calculations. Finally, to add 
perspective to the comparison of instruments described later in this thesis, a section is 
included describing the historical development of wave measuring instruments and their 
capabilities. 
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Figure 2-1: Energy Density in Sea Surface Waves as a Function of Frequency [13] 

2.1 Surface Wave Theory 

2.1.1 Wave Propagation 

Ocean surface waves behave in a predictable manner (if one overlooks their generation and 
dissipation). While the ocean surface wave height as a function of space and time is a 
stochastic process [43] (see discussion below), the kinematics of freely-propagating waves 
is governed by the dispersion relation, which relates the frequency of a surface wave, u> 
(rad/s), to its wavenumber, k (rad/cm) where k = 2ir/\ with A = wavelength (cm). For a 
linear, inviscid, irrotational flow, in the absence of currents, with both gravity and surface 
tension restoring forces, it takes the following form 

tanh(fch) (2.1) 


u, 2 = 


gk + 


rk 3 
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In deep water (kh 1) this simplifies to 


u 


gk + 


t/c 3 


( 2 . 2 ) 


where g is the gravitational constant (981 cm/sec 2 ), h is the water depth (cm), r is the 
surface tension (~ 73 dynes/cm), and p is the density of sea water ( 1 g/cm 3 ). 

When k <C 3.7 rad/cm the gravitational term in equation 2.2 dominates and the waves 
are termed ‘gravity waves’. For k 3.7 rad/cm the surface tension term dominates and 
the waves are known as ‘capillary waves’. Because the demarkation between these regions 
is not very abrupt, waves with k ~ 3.7 are usually termed ‘gravity-capillary’ waves. 

Equation 2.2 relates the wavelength of a particular wave to its frequency and, thus, re¬ 
lates the temporal characteristics of waves to their spatial characteristics. This is important 
in practice because measuring the temporal properties of surface waves is easier than mea¬ 
suring their spatial properties. It is, however, the spatial attributes of the gravity-capillary 
waves which dominate microwave electromagnetic scattering, as will be discussed in Sec¬ 
tion 2.2. Additionally, the dispersion relation is modified in the ocean environment because 
of the effects of surface currents, wave-wave interactions, and other nonlinear phenomena 
such as wave breaking. Even the linear processes, like Doppler shifting, alter the apparent 
frequencies (and wavelengths) of short waves making the transformation from the frequency 
(temporal) domain via the wavenumber (spatial) domain using the dispersion relation is of 
limited value. The significance of these difficulties will be explored further in Section 2.3 
and is one of the motivations for the work presented in this thesis. 


2.1.2 Energy Distribution 

Assuming one can model an ocean surface wave field as a random process, the power (or 
wave) spectrum of the waves is described by the spatial and temporal energy distribution. 
As described by Kinsman [37] the covariance of the surface vertical displacement is defined 
as 


p(x, r; t, t) = rj(x, t)r]{x + r,t + r) 


(2.3) 
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where t] is the sea surface vertical displacement as a function of position, x , and time, 
t, r is the separation vector between two position vectors, and r is the time separation. 
Its Fourier transform, known as the general three-dimensional power spectrum of the sea 
surface displacement, is 

4>(fc, uj\ x , t ) = (2 tt) -3 JJJ p(x, r; t, r) exp [— i(k ■ r — wr)]drdr (2.4) 

where k = (k x , k y ) is the wavenumber vector, and lo is the frequency. 

Many interesting aspects of the full three-dimensional spectrum can be visualized by 
reducing the dimensionality of the spectrum. These include the time evolution of wavenum¬ 
ber spectra for wave growth studies, the evaluation of wavenumber-frequency spectra for 
studies of free and bound waves, and the angular distribution of surface waves to study 
directionality and spreading. The reason for reducing the three-dimensional spectrum into 
a two-dimensional one in the context of the present study is that electromagnetic scattering 
theory, discussed in the next section, depends only on an instantaneous image of the sea 
surface, so the temporal variable is not needed for the scattering analysis. To this end, when 
one takes the three-dimensional spectrum from equation 2.4 and assumes a homogeneous 
and stationary wave field, 


4>(fc, cu; x, t ) 4>(fc, u) (2.5) 

and then integrates out the frequency, u, one is left with 

/ OO 

p(r ) exp [— i(k ■ r)\dr (2.6) 

-OO 

where the surface covariance becomes 

p(r) = ri(x)r)(x + r ) (2.7) 

which holds true if it is assumed that the sea surface is homogenous, implying that the 
surface covariance depends only on the relative lag between measurements (r) and not 
the absolute location of the measurement ( x ), and if it is assumed that the sea surface is 
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stationary, implying that the surface covariance does not depend on time. Both of these 
assumptions are generally justified for sea states which are not fetch limited or which have 
a steady wind field. It is with these assumptions that surface wave height array data 
collected with a buoy deployed in relatively unchanging meteorological conditions over a 
limited period of time can be considered statistically homogeneous or consistent. 

The two-dimensional spectrum from equation 2.6 describes how each wavenumber com¬ 
ponent of the surface wave field contributes to the total mean square wave energy, t] 2 where 



Numerous experiments have been conducted to measure this two-dimensional wavenum¬ 
ber spectrum. Many have measured time series wave heights at an array of locations, either 
in wave tanks or at-sea, and used the dispersion relation to transform the results into the 
wavenumber domain. For the long wave portion of the spectrum this has been very suc¬ 
cessful. Other experiments have used remote sensing instruments such as scatterometers or 
radars to try to solve the inverse problem: given a measured scattered return and a model for 
the surface scattering as a function of the two-dimensional wavenumber spectrum, among 
other things, what is the appropriate wavenumber spectrum which would reproduce the 
observed scattered return? Each experiment has helped generate a better defined model of 
the ‘true’ two-dimensional ocean wave height spectrum. 

In this model generation process a common means of modeling the wavenumber spec¬ 
trum has developed. It uses a slice of the wavenumber spectrum in the direction of the 
dominant wind waves at the spectral peak in conjunction with a spreading function which 
models the directional spreading of the energy on the surface. The result is 

*(*) = $(**> M => 9) = $(k, 9 max )D(k , 9) (2.9) 

where 9 max is the direction around which the majority of the energy in the spectrum is 
centered, and is believed to be closely aligned with the wind direction. The shape of 
the wavenumber spectrum for the short wave region in the dominant wave direction has 
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generally taken the form 


$(k,9 max ) cc k m (2.10) 

where m has taken on values from 3 to 5, depending on the experimental conditions [48, 
46, 7, 2], The proportionality constant often includes parameters such as wind speed, wind 
stress, fetch, peak spectral frequency, and/or wave age. 

The early models for the directional spreading were based on low frequency limits and 
low angular resolution using pitch-and-roll buoys such as the work of Longuet-Higgins et 
al. [40] who modeled the directional spread as a function of frequency 

J(/,8) = /lcoa 2 ‘( g ~^ al ) (2.11) 

where s was determined through a parameterized fit to the data. Later models modified the 
s-parameter, such as that of Hasselmann et al. [26], but agreement on a consistent behavior 
for s has not been achieved. Because of this inability to adequately define s, Donelan et 
al. [17] (DHH) proposed an alternate form of the spreading function. Their new model was 
based on the idea that envelope soliton groups of Stokes waves propagating around the main 
wave direction can be described by a hyperbolic secant. The energy distribution about the 
primary long wave direction should, therefore, also take the form of a hyperbolic secant, 
assuming that the wave field consists primarily of these envelope solitons and noting that 
they survive interactions with each other. DHH showed that their proposed form of the 
spectrum 


D(w, 9) = A cosh _2 (/?0) (2.12) 

where (3 is a function of the frequency ratio w/u; pea k and takes on different values in different 
w / w peafc regions, fit their capacitive wire wave gage array data more accurately than the 
previous cos 2s (#/2) form did. 

As more in situ experimental data was collected at higher frequencies and with better 
angular resolution, however, it became apparent that the DHH (3 was not universally ap- 
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plicable. Furthermore, Banner [3] disagreed with the DHH extrapolation of (3 to higher 
wavenumbers using a cutoff criterion which limited spreading at higher wavenumbers and 
found that when he used the DHH spreading function without the cutoff, he was able to fit 
both the old and newly available data more accurately with 


P={ 


2.28 (k/k p )-°- 65 

2Q—0.4+0.8393 exp[— 0.567 !n(fc/fcp)] 


for 0.97 < k/k v < 2.56 
for k/kp > 2.56 


(2.13) 


where k p is the wavenumber at the peak of the spectrum. 

Apel [1] followed in Banner’s footsteps but chose to use a Gaussian form of the spreading 
function which matches the Banner modification to the DHH spreading function to within 
experimental errors. He additionally extended its applicable range out to the gravity¬ 
capillary region based on tank tests conducted by Jahne and Riemer [34], He also used 
his resulting two-dimensional spectrum as an input to a scattering model and compared 
the results with available radar data. (Previous models had been very loose when it came 
to representing the two-dimensional spectrum either as a function of u> or k , but as data 
has become available in the higher wavenumber region, all the new models use only the 
wavenumber form for the reasons described earlier.) The spreading function Apel used is 

where 


1 


0.U + 5.0{k/k p )~ i:i 


(2.15) 


Elfouhaily et al. [19], after analyzing the Bjerkaas and Riedel spectrum [7], the Donelan 
and Pierson spectrum [18], and the Apel spectrum [1], proposed yet another form of the 
spreading function based primarily on the Apel shape. This new function, however, included 
an important aspect that the previous models did not: it required that the spreading 
function be centrosymmetric (symmetric through the origin). This property is imperative 
if the spectral model is to be used in electromagnetic scattering models (a key aspect of 
this thesis) because the electromagnetic scattering theory depends heavily on the Fourier 
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transform of the height covariance function. Because the directional spectrum is the Fourier 
transform of the two-dimensional height covariance function (as was defined previously), it 
is, by definition, a real and even function, leading to the symmetry requirements on the 
spreading function. Elfouhaily et al., thus, restructured the form of the spreading function 
so that 


D(k , 9) = -5- [1 + A (k) cos(2#)] (2.16) 

2tt 

When the spreading functions of Mitsuyasu et al., Hasselmann et al., Donelan et al., and 
Apel are modified to fit Elfouhaily’s form they all involve the hyperbolic tangent. Recog¬ 
nizing their contributions, but also their limitations, Elfouhaily et al. modeled A as 


A (k) = tanh a 0 + a p (c/cp) 2 ' 5 + a m (c m /c) 


\ 2.5 


(2.17) 


where a 0 = ln(2)/4, a p = 4, a m = 0.13u*/c m , and c m = 0.23 m/s. 

Except for Apel’s model, the previously described directional spreading models were 
based on parameterized fits to in situ experimental data. An alternate method for measuring 
the directional spectrum on the ocean surface is to use remote sensing tools and an inverse 
processing technique. This will be described in more detail in Section 2.2. Two recent 
authors have used this method and added their own modifications to Apel’s spectrum. 
Caudal and Hauser [9] argue that the spreading proposed by Banner cannot be extrapolated 
past the short gravity wavelengths. They choose, therefore, to add a modifying term, S, to 
Apel’s spectrum which functions to narrow the spectrum at higher wavenumbers. 



0.14 + 5.0(k/k p )~ 13 + 5 


(2.18) 


where his resulting form of 6 after incorporating the directional spectrum into a scattering 
model and comparing it to P and K Q -band high-frequency radar data is 


S = io (pi x2 +P2 x +ps) 


(2.19) 


where X = log 10 (fc), pi = -0.177, p 2 = 1.11, and p 3 = -2.00. 
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Finally, Romeiser et al. [55], accepted the basis of the Apel form of Donelan’s spectrum, 
but noted that it did not accurately reproduce measured cross-wind/downwind ratios mea¬ 
sured with scatterometers from L to K a -band. The model they propose takes the following 
form 


1 

2 W) 


0.14 + 0.5 


1 — exp ( — 


kuwY 
ci /. 


+5 exp 


2.5-2.6 In 


1.3 In 

V u n ) 



( 2 . 20 ) 

( 2 . 21 ) 


where ci = 400 rad/s, k n = 1 rad/m, u n = 1 m/s. 

A more comprehensive reference list of the development of spectral models over the last 
five decades is included in Appendix A. 


2.2 Electromagnetic Scattering Theory and Remote Sensing 

Remote sensing instruments, such as radars and scatterometers, generate electromagnetic 
radiation, direct it at the ocean surface, and then use the backscattered or reflected signals 
to discern information about ocean surface properties. Interpreting the reflected signals 
from the surface requires detailed analytical analysis and a thorough understanding of 
electromagnetic surface scattering theory. The theory of electromagnetic scattering from the 
ocean surface will be discussed first, and then its mathematical development and application 
using the normalized radar backscatter cross section will be explained. 

2.2.1 General Theoretical Development 

The first scattering theories were based on geometrical optics arguments which assumed 
that the reflective surface was smooth, in the sense that the radius of curvature of the sur¬ 
face was large in comparison to the impinging electromagnetic wavelengths. These theories 
were only able to accurately predict scattering at near normal incidence angles, however, 
because they modeled the sea surface as being composed only of very long waves. It was 
later found that short surface waves with wavelengths near the wavelength of the incident 
electromagnetic field contributed significantly to the backscattered return at moderate in- 
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cidence angles. To try and incorporate these small waves a new theory was presented in 
the 1950’s which described the scattering surface as being composed of small perturbations 
from a smooth, possibly undulating surface [54], This model, describing ‘slightly rough 
surfaces’, used the concept of Bragg type scattering which results from coherent or resonant 
reflection off of an ordered surface. Bragg scattering predicts that the scattered return from 
a roughened sea will be proportional to the (two-dimensional) Fourier transform of the sur¬ 
face elevation spatial covariance function at a particular wavenumber. The wavenumber is 
dependent on the electromagnetic wavelength used by the scatterometer or radar, and the 
angle of incidence of the radiated beam. Furthermore, if the sea surface is in motion, the 
scattered return will be Doppler shifted by a frequency associated with its wavelength and 
the velocity of the waves upon which the Bragg scatterers are traveling. The limitations of 
this method are the product of its assumptions of small surface displacements compared to 
the electromagnetic wavelength (inherent in all perturbation approaches) and small surface 
slopes. While a reasonable approximation, these assumptions do not always apply to real¬ 
istic sea surfaces, particularly as they become rough due to increased wind speeds. More 
refined theories were therefore formulated to remedy the discrepancies between observation 
and mathematical description. 

Research in the 1960’s showed that these limitations could be overcome if the ‘slightly 
rough surfaces’ theory was further modified by assuming that the large scale waves on the 
sea surface could be broken into small patches of tangent planes upon which the shorter 
scale scatterers propagated [4, 5, 69]. The perturbation theory, therefore, is applied to 
individual small patches that are summed to get the total scattered return, assuming that 
the short waves decorrelate very quickly relative to the large scale waves. This theory, 
known as composite surface theory or two-scale theory, has been used as a basis for most of 
the recent models of sea surface scattering [67, 18, 52, 39] despite the limitations imposed 
by its linear assumption which comes primarily from the limitation that the small patches 
must be uncorrelated. At the very large and very small incidence angles, higher order 
approximations using surface currents [56], modulation effects [25, 58], bound or parasitic 
capillary waves [53, 47], and specific wave shape features such as wave crests [36] and 
bores [64, 42] are in development but have not been applied generally in scattering analysis. 
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2.2.2 Mathematical Development 


Maxwell’s equations form the basis of electromagnetic theory. A number of texts describe 
how they can be used with the appropriate boundary conditions and assumptions to develop 
simple scattered field equations from a smooth surface for both vertically and horizontally 
polarized incident plane waves [51, 65]. These scattered field equations can then be used to 
determine the mean power received by an antenna from a surface. This is accomplished by 
integrating the scattered field equations over the antenna aperture weighted by the antenna 
beam pattern, and then normalizing appropriately. Using the slightly rough surface theory 
assumptions (and assuming the sea surface is homogeneous) it can then be shown that 
the non-dimensional normalized radar cross section of a slightly rough saltwater surface at 
moderate incidence angles can be expressed as [50] 

a 0 (<M) = Stt^G^TOM) (2.22) 

where ko is the electromagnetic wavenumber in units of radians/length and G(4>) is a non- 
dimensional polarization dependent geometric reflection coefficient incorporating the dielec¬ 
tric constant such that 


Gh{4>) = 


cos 4 ((f)) 

[1 + (l/Ve)cos(<£)] 4 


and 


(2.23) 


= cos-M(l + si^W 
(cos(0) + (1/v^)] 4 


(2.24) 


where (f) is the incidence angle in radians, 9 is the azimuthal look direction in radians, e is 
the non-dimensional relative dielectric constant of sea water [57], and finally 


T(0, 9) = $(2 k 0 sin(</>), 9) + $(2 k 0 sin(<£), 9 + n) (2.25) 

which has units of length 2 /(radian/length) 2 . These equations show that the normalized 
radar cross section from an air-sea interface using the Bragg backscattering model is directly 
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proportional to the wavenumber spectrum evaluated at the Bragg resonant wavenumber, 
k B = 2kosm((f)). 

As was noted in Section 2 . 2 . 1 , when the small slope assumptions are violated, as they are 
when the small scale roughness ‘rides’ on longer swell, the pure Bragg model can no longer 
be used. A modified Bragg model called the two-scale or composite surface theory is used 
instead. This model takes into account the long wave slope through a probability function 
as well as altering the apparent incidence and azimuthal look directions at which the terms 
are evaluated. Mathematically it can be derived from the simple Bragg model by taking a 
few simple steps. To begin with, the plane of incidence must be tilted to incorporate the 
long wave slope. Let the long wave slope in the plane of incidence be 7 l v and out of the 
plane of incidence be 7 Bx . The new local incidence angle then becomes 

</>' = <!>- 7 by (2.26) 

Furthermore, tilting out-of-plane is equivalent to rotating the tilted surface by a small 
angle O’ where for small angles cos O' m 1 and sin O' « ^ Lx / sin<// [69]. Now, for a vertically 
polarized incidence field 


E — E sin O'a + E cos 9'b (2.27) 

where a is a unit vector lying in the tilted surface and perpendicular to the incident elec¬ 
tromagnetic field vector and b is a unit vector lying in the local plane of incidence and 
perpendicular to the incident electromagnetic field vector. The resulting incident vertically 
polarized electromagnetic field is, thus, a combination of both vertical and horizontal po¬ 
larizations. The same can be shown for the horizontal incident field. Propagating these 
modifications through the slightly rough surface equations generates the normalized radar 
backscatter cross section for the composite surface theory 

<7 0 ((f>, 6) = 87 Tko J G(<f>, 7 Lx, r rLy)T(4>', e')P{~ 1 Lx , -lLy)d^LxdlLy ( 2 . 28 ) 

where P{jLx, iLy) is the probability distribution of the long wave slopes, and it is generally 
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assumed that T(<(>', 9 1 ) can be approximated by T(<p',6). Additionally, the G(0, 
function is defined as 


and 


Gll(<t>,lLx,7Ly) = 



(2.29) 


Gv (0, iLx, 7Lt/) 



Gvtf) 


(2.30) 


since for most incidence angles Gv » G H - Prom these last two equations it can be seen 
that the tilt modulation induced by the long wave slopes affects the horizontally polarized 
return much more than the vertically polarized return. For further detail and explanation 
of this derivation see reference [51]. 

Plant [50] went a step further with this derivation by using the Cox and Munk [14] 
Gaussian approximation of the slope probability distribution and evaluated the integral in 
equation 2.28. The results were simplified by Lyzenga [41] assuming a long wave zero mean 
slope and presented in the following form 


oo(<A 0) = cr o (0, 9) 


1 + y ('lly) + y<7L> 


(2.31) 


where 


Ul _ J_ (, ? G"(0)r'(0, b ) r"(0, 9 ) 

a 0 \d^ ) G( 0 ) + + 


(2.32) 


and 


J _ (^£^1 - 2 ( 2 v ^~ 1 ) 

ctq ^ da 2 J sin 2 ( 0 ) 


(2.33) 


where (q£ y ) is the in-plane long wave slope variance, ( 7 2 X ) is the out-of-plane long wave slope 
variance, R = Gv/Gh for horizontal polarization, R = Gh/Gv for vertical polarization, 
and a = tan -1 ( 7 ^) which is the tilt angle perpendicular to the plane of incidence [41]. 
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Equation 2.31 is an approximation to second order in the tilting angles which can be used 
to estimate of the influence of long wave slope variations on the normalized radar backscatter 
cross section. 

Other non-linear and non-Bragg mechanisms which have been postulated to affect the 
radar backscatter cross section will not be addressed in this thesis, but a detailed refer¬ 
ence list of the development of scattering models over the last five decades is included in 
Appendix B. 

2.3 Wave Measuring Instruments 

Many instruments have been developed to measure sea surface spectral characteristics to test 
and provide data to verify scattering theories. The variety of methods used includes aerial 
photogrammetry based on sun glint [14], stereo photography based on surface slope [28], 
resistive wire wave staff measurements [38], capacitive wire wave staff measurements of 
surface height [45, 17], radar and scatterometry measurements [59, 68, 21], single point laser 
slope measurements of surface height [63, 29], scanning laser slope measurements [70, 32, 8], 
and surface imaging devices using video processing [33, 71]. A more detailed review of 
the development of wave measuring instruments over the past five decades is included in 
Appendix C. 

Each of the techniques described above has its drawbacks. The wire wave gage tools 
suffer from problems such as disturbing the wave field itself, limited resolution of wave 
height or wave slope, and difficulties accommodating the many orders of magnitude that 
wave heights span (i.e. long waves with amplitudes of meters and capillary waves with 
amplitudes of millimeters). The laser devices can generally only measure slope (although 
some single point devices can be made with laser ranging capabilities), and many have 
limitations on the time of day they can be deployed due to the interference of sunlight with 
optical sensors. Additionally, the laser systems do not deal well with the change in focal 
length caused by water surface height variation over the sensor area. The remote techniques 
(photogrammetry, stereo photography, radar and scatterometer methods, surface imaging) 
suffer from difficulties in absolute calibration, and, in the case of scatterometry, lack of full 
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understanding of electromagnetic scattering from the ocean surface. All extant sensors also 
have difficulties or outright failures when measuring cresting or breaking waves. 

The resistive circular wire wave gage described in Chapter 3 is not exempt from the limi¬ 
tations listed above. It is, however, different from the instruments referenced in Appendix C 
in that it combines the simplicity of wire wave gages with the higher spatial and directional 
resolution of recent laser scanning slope gages. This directional resolution combined with 
the inadequate understanding and experimental data to describe the spreading function 
described in Section 2.2 makes the comparison of coincident data collected with both the 
circular wire wave gage and the X-band scatterometer instruments a unique opportunity to 
improve the current state of knowledge in this area. 
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Chapter 3 


Circular Wire Wave Gage Buoy 


The circular wire wave gage buoy and its data processing algorithms were designed to gen¬ 
erate an estimate of the two dimensional sea surface spectrum. Details of the prototype 
design and testing of the circular wire wave gage can be found in the author’s Masters 
thesis [61]. The prototype gage was designed as a proof-of-concept instrument and suc¬ 
cessfully achieved this goal. It was not, however, constructed with the attributes necessary 
for its proposed use in the scientific study described in this document. The study involved 
experiments conducted at sea which have formed the fundamental contribution of this re¬ 
search to the fields of oceanographic remote sensing and surface wave studies. The primary 
deficiencies of the prototype instrument included the fragile nature of its structural design, 
the need for an external deployment platform, inconsistent electrical performance, the lack 
of an external data acquisition system, and finally, the absence of an internal power source. 

This chapter details the development of the circular wire wave gage buoy beyond its 
prototype phase. It describes the system layout, the instrument’s basic components, how 
it is operated, how the data is processed, as well as the major modifications made which 
facilitated its use in the experiments described in Chapter 4. 


29 




Figure 3-1: Components of the Circular Wire Wave Gage Buoy 

3.1 Overview/Introduction 

3.1.1 System Layout 

The system, including both the instrument and data processing algorithms, is comprised 
of five major components: the filtering package, the measuring device, the data acquisition 
system, the buoy structure and housings, and the data processing package. The relation¬ 
ship between these components is schematically represented in Figure 3-1. Each of the 
components is described in detail in this chapter. 

3.1.2 Fundamental Operating Principle 

The key component of the wave wire buoy is its measuring device: the 64-element wave wire 
array. Each of the 34-gauge nickel chromium (nichrome) wave wires has a resistance of 17 
ohms/foot and acts as a voltage divider as shown in Figure 3-2. The relationship between 
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Figure 3-2: Voltage Divider Relationship for Wire Wave Gage 


the input and output voltages of a voltage divider is mathematically represented by 


V, 


R 


out 


■gage 


R + R, 


-V ir 


■gage 


(3.1) 


where R is a fixed resistance internal to the wave gage circuitry and R gag e is the resistance of 
the wave wire. If R is significantly greater than R gag e, this relationship can be approximated 
by 


V, 


R, 


out 


■gage 

R 


v ir , 


(3.2) 


Because R is fixed, the output voltage will be nearly linearly related to the input voltage 
through R gage ■ The seawater, due to its high conductivity, effectively acts as a short circuit 
between the point where each resistive wire pierces the water surface, and the aluminum 
center post of the gage which is connected to the electrical ground of the electronics package. 
As the level of the seawater on each of the wires changes, the effective resistance of each 
of the wires, R ga ge, changes, varying the voltage measured across the voltage divider. The 
voltage measured across each wire is, thus, linearly related to the sea-surface elevation at 
that point. This linear relationship allows the filtered output signal from the wave gage 
to be used to reconstruct the input surface elevations. These surface elevations are then 
used to generate an estimate of the directional sea surface spectrum which is the ultimate 
product of the instrument and its processing algorithms. 
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3.2 Buoy Structure and Housings 


The circular wire wave gage mechanical de¬ 
sign was based around the original circular 
wire layout arranged by former MIT gradu¬ 
ate student, David Chen which arranged the 
64 wire wave gages in a 20 centimeter diam¬ 
eter circle [12]. This layout facilitates the 
data processing and is fairly straightforward 
to implement. The remainder of the buoy 
design was chosen primarily to support and 
locate the wave wires, and additionally to 
minimize the flow disturbance around them 
while maintaining a stable long wave follow¬ 
ing platform. A number of additional con¬ 
straints on the design were also incorporated. 
These included low cost, due to limited fund¬ 
ing, and ease of deployment, assembly, and 
disassembly, because of intended single or 
two-person use and testing. The resulting 
buoy design is shown in Figure 3-3. 

The water-tight central housing contains 
the pre-processing analog electronics, com¬ 
pass, and temperature probe. The structure 
also functions as a support and locator for 
the 64 nichrome wires. It is constructed en¬ 
tirely of aluminum and is shown in Figure 3- 
4. 

The three canisters which form the base 
corners of the pyramid of the buoy contain 
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Figure 3-3: Top view and partial side view of 
the assembled circular wire wave gage buoy. 
Note that the wave wires (not shown) sur¬ 
round the vertical post of the central filter 
housing. 
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Figure 3-4: Central Analog Electronics Housing 
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the data collection system and the batter¬ 
ies. Up to four 12-volt rechargeable lead- 
acid batteries are secured in each of the two 
battery mounting brackets and loaded into 
two of the canisters which are constructed 
of PVC and are shown in Figure 3-5. This 
allows the flexibility of using up to eight bat¬ 
teries depending on the particular deploy¬ 
ment specifications (e.g. duration of deploy¬ 
ment, weight considerations, etc.). A vent 
port in the lid of each battery canister al¬ 
lows for any out-gassed fumes from the bat¬ 
teries to be discharged between buoy deploy¬ 
ments and during battery recharging. The 
third canister houses the PC104 form factor 
data acquisition system and is constructed 
of aluminum. The size and weight of each of 
the three housing with its contents are nearly 
equal (when using four batteries in each of 
the battery canisters), reducing the need for 
additional ballasting. The mounting brack¬ 
ets used to secure the PVC and data acquisi¬ 
tion canisters to the buoy are shown in Fig¬ 
ure 3-6. Hose clamps are run through the 
vertical oval holes and around the canisters 
to keep them from moving after assembly. 
Figure 3-7 displays the data acquisition can¬ 
ister endcap and describes its ports. 



Figure 3-5: PVC Canisters for Batteries 



.“N 






0 0 ! 




| 0 0 




j, 

1 . 



Figure 3-6: Canister Mounting Brackets 
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Figure 3-7: Data Acquisition Canister Endcap. Proceeding clockwise from the top center 
connector: Serial Port 4-pin connector, Unused 7-pin connector, Ethernet 10Base2 connec¬ 
tor, Battery 4-pin connector (one 2-pin male connector comes from each battery housing), 
and in the center, Analog Electronics Output 76-pin connector. The corresponding dummy 
plugs and watertight covers are shown next to each connector. 

The superstructure is constructed of aluminum and fiberglass angle. The angle is bolted 
to the canisters and the central housing during assembly. While a fixed/welded structure 
would have been more rigid, the bolted construction was chosen because it would have 
been awkward to transport and work on the assembled instrument due to its large size 
(2.1 meters on a side). The resulting modular nature of the super-structure allows easy 
assembly, disassembly, and transport. 

When deployed, the buoy floats with about half the length of the nichrome wave wires 
exposed above the mean water level. It follows the long waves (waves with wavelengths 
greater than 5 or 6 meters) very smoothly, and has almost no vertical motion when only 
high frequency (very short) waves are present. In choppy seas with 2 to 5 meter intermediate 
steep waves, however, it has difficulty keeping the floatation canisters from dipping below 
the water surface because of its light weight and broad footprint. For these reasons it can 
be deployed both in the open ocean and in coastal waters, but should not be deployed in 
seas where short steep waves or whitecapping is occuring. 
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3.3 Analog Electronics (Filtering Package) 


Many modifications were made to the original pre-processing analog electronics of the instru¬ 
ment. A number of these were based on the prototype performance such as improvements 
in noise reduction, changes to the analog filter shape, and reduction of analog system power 
draw. They were originally secondary to the fundamental matter of simply building a buoy 
which functioned mechanically, and building the new stand-alone data acquisition system. 
Once the new system was assembled it worked well in the lab, but nearly half the channels 
were severely attenuated when it was deployed in the ocean. Because of the difficulty in 
establishing the cause of the attenuation, some of the originally proposed modifications to 
the analog electronics and numerous new changes were implemented to solve this problem. 

3.3.1 Layout 

The analog electronics are distributed on five PC boards which are stacked horizontally on 
top of one another. Figure 3-8 shows this configuration. Four of the boards are identical 
and each contains sixteen identical channels of analog filtering. The fifth board contains the 
power conditioning, the calibration switch, the flux gate compass, the temperature probe, 
and acts as the interface between the filtering boards and the data acquisition package, 
and between the filtering boards and the wire wave gages. The electrical diagram for this 
interface board is shown in Appendix D. 

3.3.2 Filter and Modifications 

The generation and filtering of each of the 64 channel voltage outputs is accomplished in 
several steps which are illustrated in the electrical circuit diagram shown in Figure 3-9. 
First, a 10 KHz oscillating voltage is generated. This carrier signal is common to all the 
channels. The signal is then passed through a gain adjust, and sent to each of the nichrome 
wave wires. The water height on each wire amplitude modulates the high frequency carrier. 
The output of the wave wire voltage divider circuit discussed in Section 3.1.2 is then passed 
through a gain adjust and then through an RMS-to-DC converter which filters out the 
carrier frequency and leaves the envelope (the amplitude modulated signal). The signal is 
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Figure 3-8: Image of Analog Electronics Boards Stacked for Installation in Buoy. The large 
distance between the lowest board and the four filtering boards is to accommodate the 
compass. The upper copper shield board is not shown. 
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next sent through a trim offset and gain adjust, and then through a bandpass filter. The 
gain of the signal is again adjusted, and it is sent through a noise rejection low pass filter 
before finally being sent to the analog-to-digital converters in the data acquisition system. 

Modifications made to the original analog electronics are readily apparent when compar¬ 
ing Figure 3-9 with Figure 3.5 from the author’s Masters thesis [61]. The key differences are 
the addition of capacitors across the feedback loops of the filters, new operational amplifiers, 
and the modification of the filter cutoff frequency and shape accomplished by modifying 
the resistor and capacitor components of the filters. The manual calibration switch has also 
been changed to a relay controlled switch, allowing it to be computer controlled. Details 
about the calibration switch and oscillator input can be found in Appendices D and E. 

The system function of the filter electronics, representing the ratio between the Laplace 
transform of the response of the system to the transform of the source [60], can be mathe¬ 
matically estimated using the nominal values of the circuit elements. Table 3.1 shows the 
circuit diagrams and resulting mathematical forms of the system functions for the individual 
steps in the filtering process. Combining all the individual filter system functions by multi¬ 
plying them together gives the system function for the entire bandpass filter. Substituting 
in the appropriate element values and evaluating along the frequency axis (s = jw) results 
in the sinusoidal steady state frequency response of the analog filter for each channel. A 
plot of this appears in Figure 3-10. The filter shape was chosen primarily to whiten the 
input energy in the low frequency range. This is necessary because as can be seen from 
Figure 2-1, the region of the spectrum from 0.3 Hz out to about 2 Hz includes nearly four 
orders of magnitude in energy. In order to resolve this vast range with a data acquisition 
system having limited resolution it was necessary to whiten the input. It is also impor¬ 
tant to note that given the Nyquist frequency for the data sampling of 50 Hz, it is evident 
that the analog filter is adequately designed to prevent any significant aliasing from high 
frequency input or noise. 

3.3.3 Noise Reduction and Signal Conditioning 

While the noise in the prototype analog electronics was acceptable, since a new housing and 
power control board were being developed, some effort was spent trying to reduce the noise 
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Figure 3-9: Single Channel Front End Analog Electronics 
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Name Circuit Diagram 
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Table 3.1: Circuit Diagrams and System Functions for Analog Filter 
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Figure 3-10: Frequency Response for Analog Filtering of Each Channel. The short dashed 
lines represent the approximate limits of the digital bandpass filter used in the post process¬ 
ing. The long dashed line represents the Nyquist frequency of the data acquisition system. 
The x-axis is shown in log scale for the lower two plots. These functions are very repro¬ 
ducible from channel to channel. 
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even further. To this end, the following modifications were made to the analog system. 
Some of these modifications were made to help solve the attenuation problem mentioned 
earlier. 

1. Capacitors were added between the plus and minus voltage supplies and ground on 
all the integrated circuit (IC) chips. 

2. New amplifier ICs were used which reduced the quiescent current draw by each am¬ 
plifier by more than a factor of ten. 

3. Plain PC boards with a thin sheet of copper which are the same physical dimensions 
as the electronic boards were mounted on either end of the electronics stack and 
grounded to the housing. 

4. A ground wire was added to each of the four electronics boards connecting it directly 
to the housing. 

5. The cable running from the analog electronics housing to the data acquisition system 
was fully shielded. 

6. A temperature probe was placed in the analog electronics housing to monitor the 
environmental temperature of the capacitors so that any temperature changes which 
might affect the value of the capacitors (which are temperature sensitive) can be taken 
into account during calibration. 

3.4 Data Acquisition Hardware 

The data acquisition system was developed to sample and store a set (or scan) of the 64 
analog outputs at 100 Hz (i.e. the beginning of each scan is triggered every 0.01 seconds, 
but the collection of the 64 analog voltages occurs much faster than that). It also collects 
data from the on-board compass, the batteries, and the temperature probe. It was designed, 
built, and programmed by Kenneth Prada, at the Woods Hole Oceanographic Institution. 
The specifications and performance of the system are detailed in Table 3.2. A complete list 
of the design criteria and constraints is presented in Appendix F, and a schematic of the 
data acquisition wiring diagram can be found in Appendix G. 
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Stack configuration 

PC/104 

Processor 

CoreModule 486-11 50 MHz 

A/D Boards (4 each) 

16 channel RTD DM5406 

Network board 

Ampro MiniModule/Ethernet 

Hard drive 

Western Digital 540 MB, Model WDAL2540 

Batteries (8 each) 

Yuasa NPH5-12, 12V 5.0Ah rechargeable 

Data acquisition speed (per scan) 

100 Hz 

Burst rate for sampling purposes 

100 kHz 


Table 3.2: Data Acquisition Basic Specifications 


3.5 Power 

The wave gage is powered by up to eight 12 volt rechargeable batteries. These batteries are 
mounted in two wooden storage racks which fit snuggly into the PVC housings to prevent 
movement during deployment. The four batteries in each housing are wired in parallel 
with appropriate diodes to allow for recharging and prevent them from discharging each 
other. The wiring diagram for this arrangement is shown in Figure 3-11. The batteries 
are typically charged between each deployment. To do this, the vent port in each battery 
housing is opened and the positive terminal of the charging circuit, which is attached to 
the vent plug, and the ground connector on the battery output cable are connected to a 
battery charger. The maximum current which should be used is 0.5 amps per battery, and 
the maximum voltage charge should be set at 15.7 volts. This will charge the batteries to 
13.8 volts, their maximum rated charge, and will allow for losses in the system due to the 
diodes and load resistors. After a full day’s deployment the batteries may need to charge 
for up to ten hours. 

The entire system (digital and analog) draws approximately 19 watts when running. 
With eight batteries fully charged and a safety factor of two, the instrument can collect 
over 12 hours of data. In order to conserve power during a deployment, a Watchdog timer 
is used which shuts down the power to all the systems between data sets, and then restarts 
them at a user specified interval, or until a user specified number of data files has been 
collected (see Section 3.6.1). This power save mode allows the instrument to be deployed 
for longer durations than it would be limited to if it were collecting data continuously. It 
was found, however, that in productive coastal waters the wires became coated with film 
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Figure 3-11: Battery Wiring and Charging Diagram 

from the water surface, which affected their ability to accurately measure wave height after 
only about four to five hours, so the deployment duration is not generally controlled by 
power limitations. 

3.6 Operation 

The wave gage buoy is controlled by two programs. The primary program, wave.c, is the 
data acquisition code which collects a single data set including signal outputs from the 64 
wave wires, the temperature probe, the compass, and the battery voltage. A batch file, 
run_wave.bat, contains the call to wave.c and its command line options. The secondary 
program, setup, c, sets the start time for the first call of run_wave, the interval between 
the beginning of each call to run.wave, and the total number of data sets to be collected. 
Copies of wave.c and setup.c can be found in Appendix H. Once all the data files have 
been collected, the system shuts down and does not ‘wake up’ again. Once the data has 
been collected it can be downloaded through a local network which can be accessed via a 
10Base2 port on the data acquisition housing. 

3.6.1 Startup and Data Download 

The data acquisition (daq) system is turned on by plugging the power cords from the two 
battery housings into their connector on the daq housing endcap. The computer system 
operates in two basic modes: data collection mode and serial console mode. In the data 
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collection mode, which is the default, the system turns on, collects data based on the user 
specified startup configuration settings, and then goes to sleep (shuts down all power). In the 
serial console mode the user must be connected to the buoy through a serial cable and have 
a serial console program running on the externally connected computer. When the system 
turns on, the user is given three choices for operation, the ‘operate’ state (data collection 
mode), the ‘normal’ state (dummy terminal access to the wave gage buoy computer, i.e. 
C:\ prompt), or the ‘PC-Anywhere’ state, which runs the setup programs on the buoy to 
allow it to be accessed using PC-Anywhere as an alternate means of transferring data to 
and from the buoy. Data can be downloaded via a 10Base-2 network connection when one 
is in the serial console mode running in the ‘normal’ state. 

Startup 

The user must do two things to set the buoy up to collect data autonomously. Both must be 
done in the serial console mode. First, the length of each of the data sets must entered into 
the run_wave.bat program by setting the number of blocks to collect (see Section 3.6.2). 
This batch file contains only one line: 

wave -b ??? 

where the questions marks are replaced with the appropriate number of blocks to collect. 

The second thing the user must do is run the setup program. The setup configuration 
program is accessed by running the setup. exe program from the C: \ prompt. A description 
of the program output and expected user input is detailed in Table 3.3. Once the settings 
have been stored, the system powers down and waits for its initial start time to arrive. If 
there is a problem with the startup configuration, it can be reset by unplugging the power 
from the daq system for at least fifteen seconds. This resets the Watchdog timer. Upon 
reconnection of power the system will run in its default data collection mode unless a serial 
console is running and a user selects the normal state as described above. The user will 
then need to rerun the setup program. 

The serial console mode is established by first, installing the serial cable between the 
buoy and a remote computer, second, running a networking software package on the remote 
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Preamble 


WARNING!!!! 

This program performs a system setup 
Then it shuts the system down with alarm set! 

Press <RETURN> to continue or <ESC> to quit.... 


Current time 
User input 


Confirmation 
Shut down 


Current Time: hh.mm.ss 


Enter Start Hour.? < hi > 

Enter Start Minute.? < ml > 


Enter Number of Data Sets..? 
Next Start Time: < hi >:< ml > 
Shutdown in 2 seconds. 


Table 3.3: Setup Configuration and User Prompted Information 


computer which allows serial connections, and third, powering up the daq system. The serial 
cable is a 4-conductor cable which plugs into a water-tight connector on the daq housing 
endcap and the serial port on the remote computer. The networking software package must 
be configured as described in Table 3.4. Hyper Terminal made by Hilgraeve works well as 
a simple serial console interface program for this purpose. It should be noted that only 
programs on the daq system which can be spooled through the serial port can be run on 
the remote computer when operating in serial mode. For instance, the text editor edit. com 
cannot be run in serial console mode, but the text editor me. exe (Micro Emacs) can be. 


Data Download 

The data stored on the wave gage buoy hard drive can be easily downloaded to another 
computer without having to open the daq aluminum housing. This is accomplished by using 
the following procedure: 

1. Connect a 50 f l BNC (10Base-2) cable between the data acquisition endcap and the 


Baud Rate 
Data Bits 
Parity 
Stop Bit 
Flow Control 

9600 

8 

None 

1 

None 

Terminal Emulation 

ANSI 


Table 3.4: Serial Console Configuration 
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downloading computer network port with 50 tt terminations at each end 


2. Start a serial console mode as described above 

3. Change directories to C:\data on the buoy computer 

4. Run the startnet batch file on the buoy computer 

5. Run the ftpsrv program on the buoy computer 

6. ftp to 128.128.93.22 (the buoy IP address) from the downloading computer 1 

7. Set the transfer mode to binary 

8. Download the appropriate data files and exit ftp when finished 

9. Exit the ftpsrv program on the buoy by typing ‘q’ at the buoy computer prompt 

10. Restart the buoy computer before using it to collect data again (this clears the network 
settings from the memory, to enable subsequent data acquisition) 

3.6.2 Basic Commands, Calibration, and File Formats 
Basic Commands 

The data collection mode of operation is entirely controlled through the wave. c program. 
This program has various command line options which are detailed in Table 3.5. An ap¬ 
propriate use of the command line options for testing purposes would be 

wave -b 100 -s 

which starts collecting 100 * 5.12 = 512 seconds of data (51200 data samples per channel) 
immediately and does not shut down at the end of the data collection. When planning to 
use the power saving mode which shuts the system down and then wakes it up at a later 
time, the run_wave. bat program along with the setup. c programs must be used instead of 
the basic command line call (as described in the Startup section). It should be noted that 

’The downloading computer must be configured with a 128.128.*.* IP address or the buoy computer will 
not recognize it 
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-g <gain> 

-i <interval> 
-b <nblks> 

-q 

-s 

-h <hour> 

-m <minute> 

-c 

-? 


Set A/D board gain 

Set time between beginning of data sets (in min.) 
Set number of 5.12 sec. data blocks to collect 
Sets quiet mode (turns off debug screen displays) 
Inhibits shutdown of daq after data collection 
Set start hour to (using 24 hour clock) 

Set start minute to 

Enable calibration step after first block of data 
Prints a help menu for the command line options 


(default = 1) 
(default = 15) 
(default = 32) 


(no default) 
(no default) 


Table 3.5: Command Line Options for wave. c 


while the command line options appear to accomplish the same purposes as the setup. exe 
program, they should generally not be used except when testing the instrument without 
using the power saving shut down mode. The reason for this is that each time the buoy 
‘wakes up’ from the power save mode to collect a subsequent data set, it references the 
run.wave.bat batch file, NOT any original wave.exe program call that may have been 
entered at the C: \ prompt. 


Obtaining Calibration Files 

The calibration files are used in the post processing of the buoy data and should be collected 
both before and after a deployment while the buoy is out of the water. Their purpose is 
to provide information about the noise in the system, to verify that each of the channels 
is functioning properly, and to determine the actual transfer function for each channel. 
The instrument is calibrated by electronically triggering a relay, which steps the amplitude 
of the oscillator voltage, and then stores the response of the system to this step. The 
electronic triggering of the step is controlled in the wave. c program when the -c command 
line option is used. To improve the statistical validity of the calibration results, eight files 
are collected during a standard calibration: four calibration files and four non-calibration 
(noise) files. This process is automated with a batch file that repeatedly calls the wave. c 
program with the appropriate command line options. This batch file, calib.bat, is located 
in the C:\calib directory. The eight files generated have the same format and file names 
as the typical wave, c output files, but they are only 2 blocks (10.24 seconds) long, and they 
are alternately a noise file, then a calibration file, then a noise file, and so on. 
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File Description and Formats 

The output data files are stored on disk with month,day,hour,minute.* file names 

mmddhhmm. wav - wave wire data from all 64 channels (stored in binary) 
mmddhhmm. tmp - temperature data (stored in ASCII) 
mmddhhmm. kvh - compass data (stored in ASCII) 

The *. wav file has a 64-byte header and the data is double-multiplexed as two-byte 
signed integers with 8 data values per line. The data is stored in 5.12 second blocks, is 
indexed by board, and is also indexed by time step; the first data point (time tl) for the first 
16 channels (board 1) is stored, then the second data point (time t2) for the first 16 channels 
is stored, and so on until 512 data points for each of the 16 channels has been stored. Each 
of the remaining three boards (channels 17-32, channels 33-48, and channels 49-64) is stored 
the in the same fashion. Once the first 512 samples of each of the 64 channels have been 
stored, then the cycle repeats for the second block of 512 samples. Figure 3-12 illustrates 
this organization. The wav2asci. exe program converts the binary <run numeric> .wav file 
into a <run numeric>. asc ASCII file and is described in Section 3.6.3. 

The *. tmp file has a 64-byte header and each data value, which has units of volts, must 
be multiplied by a calibration factor of 10° C/volt to get an output in degrees Celsius. The 
temperature probe value is not collected at the 100 Hz data collection rate of the wave 
wire data because it does not vary over 0.01 second time intervals, rather it varies over 
five to ten seconds. It is instead collected once every 5.12 seconds (after each ‘block’ of 
data). It’s primary use is to monitor the temperature sensitive capacitors in the central 
analog electronics housing. If any significant temperature shifts (more than 5°) occur, the 
capacitor values in the calibration calculations can be modified. 

The *. kvh file has a 64-byte header and each data value is stored in KVH format as 
described in Table 3.6. The compass data is not collected at 100 Hz, because of its rather 
coarse resolution (5°) and because it tends to vary only every few seconds. It is, therefore, 
collected once every second. If the seas are fairly choppy, the compass output can be very 
noisy because the gimbal on which it is mounted cannot compensate for the full motion of 
the buoy. In these cases, tests in the lab showed that computing a running window average 
over 9 seconds was a fairly accurate way to reproduce the approximate compass heading of 
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board 4 board 3 board 2 board 1 


chl(tl) ch2(tl) ch3(tl) ch4(tl) ch5(tl) ch6(tl) ch7(tl) ch8(tl) 
ch9(tl) chlO(tl) chi l(tl) chl2(tl) chl3(tl) chl4(tl) chl5(tl) chl6(tl) 

chl(t2). 


ch 1 (t512) . 


chl7(tl). 

chl7(t512). 


ch33(tl). 

ch33(t512). 


ch49(tl). 

ch49(t512). 


<N 


p o 

-c 


-C 


. chi6(t2) 


. chl6(t512) 


.ch32(tl) 

. ch32(t512) 


.ch48(tl) 

. ch48(t512) 


.ch61(tl) 

. ch61(t512) 



Figure 3-12: Data Format for *.wav File 


BLOCK 1 


BLOCK 2 


BLOCK 3 
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EXAMPLE DATA OUTPUT STRING: $abc.d,e,ff*CHECKSUM<cr> 


abc.d compass heading to tenth of a degree 
e output mode (D if in degrees, M if in mils) 

f f status of magnetic surroundings (OL=overload, OK=data valid) 

* indicates beginning of 2-byte checksum of message 

CHECKSUM hexadecimal checksum of all characters between $ and * 

Table 3.6: KVH C100 Compass Data Format 

the instrument. 


3.6.3 Diagnostic Programs 

Numerous diagnostic programs were written to allow trouble shooting and testing of the 

wave gage buoy operation and data acquisition. They are located in the c:\wave directory 

on the wave gage buoy computer. The key programs are described below. 

test_*. c Samples the 16 channels of board *, displays them, then recaptures them each 
time <enter> is hit. <ESC> exits. 

min_max.c Using ‘min_max cboard (l-4)> <channel (0-15)>’ the user selects a board and 
channel and program returns a hexadecimal and decimal display of the minimum and 
maximum voltage output from that channel over a 10 second interval. 

chantest. c Using ‘chantest cboard (l-4)> Cchannel (0-15)>’ the user selects a board and 
channel. The program returns a continuous screen display of the voltage output from 
that channel in both hexadecimal and decimal units. <ESC> exits. 

test-64, c Outputs a running screen display of all 64 wave wire voltage outputs. <ESC> 
exits. 

wav2asci.c Converts binary mmddhhmm.wav file into an ASCII file, mmd dhhmm .asc, with 
the format described in Figure 3.6.3 

battery. c Returns a continuous display of the current battery voltage. <ESC> or <RETURN> 
exits. 
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scan 1 


scan 2 
scan 3 


Figure 3-13: Output File Format of wav2asci.exe Program 

3.7 Data Processing 

3.7.1 Process Description 

One of the primary difficulties in obtaining directional spectra of the sea surface is the 
physical difficulty in obtaining a spatially two-dimensional measurement of the sea surface 
elevation. The input signal must also be measured over time at every sampling point to 
allow temporal averaging if the spatial dimensions of the region are not expansive. One of 
the features of the analysis used to process the circular wire wave gage data is that evenly 
spaced data sampling points are not required. Instead of using an evenly spaced grid, a 
more modest spatial sampling is used. The following analysis describes the method used to 
obtain the directional spectrum of the sea surface given that the spatial two-dimensional 
vertical surface covariance can be calculated at some known spatial lags. 


chl(tl) ch2(tl) ch3(tl) ch4(tl) ch5(tl) ch6(tl) ch7(tl) ch8(tl) 
ch9(tl) chlO(tl) chll(tl) chl2(tl) chl3(tl) chl4(tl) chl5(tl) chl6(tl) 
chl7(tl) chl8(tl) chl9(tl) ch20(tl) ch21(tl) ch22(tl) ch23(tl) ch24(tl) 
ch25(tl) ch26(tl) ch27(tl) ch28(tl) ch29(tl) ch30(tl) ch31(tl) ch32(tl) 
ch33(tl) ch34(tl) ch35(tl) ch36(tl) ch37(tl) ch38(tl) ch39(tl) ch40(tl) 
ch41(tl) ch42(tl) ch43(tl) ch44(tl) ch45(tl) ch46(tl) ch47(tl) ch48(tl) 
ch49(tl) ch50(tl) ch51(tl) ch52(tl) ch53(tl) ch54(tl) ch55(tl) ch56(tl) 
ch57(tl) ch58(tl) ch59(tl) ch60(tl) ch61(tl) ch62(tl) ch63(tl) ch64(tl) 


chl(t2) 


ch64(t2) 


chl(t3) 


ch64(t3) 
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The Vertical Surface Covariance 


Let the sea surface height covariance be defined as in Equation 2.3. The covariance function 
is assumed to occupy a region in covariance space bounded by 

x 2 + y 2 < R 2 (3.3) 

To simplify the implementation of the spectral estimation method to be used, a rectangular 
(or square) region in covariance space is adopted as opposed to a circular one. To accomplish 
this, additional zero-valued data points have been added to data set in the corners of a box 
bounded by 


|x| < L x and |j/| < L y (3.4) 

where r = (x, y ), L x and L y are the horizontal and vertical distances spanned by the spatial 
differences in the sample locations, and R is the radius of this occupied region (R. = L x = L y 
for a square region). A smoothing window is used such that p(r\t), the height covariance, 
smoothly approaches zero at the boundaries of this region and outside the region it is zero. 
A smoothing window is used because the covariance function is only known within this finite 
circular region, and in a similar fashion that one-dimensional smoothing windows are used 
to reduce side lobe contamination in the Fourier transform of a finite time series records, 
a two-dimensional window is needed to reduce the side lobe contamination in the resulting 
transform of the covariance data. This is accomplished with an energy conserving window, 
W, such that given 

y/x 2 + y 2 
W ~ R 

then 


W = { 


2/7r i 


arccos(w) — w\/\ — vf z 


if w < 1.0 
otherwise 


(3.5) 
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so that 


Pwin{r-,t ) = p(r;t) ■ W 


(3.6) 


and 



Wdw = 1 


(3.7) 


Let the windowed covariance function be represented by an infinite sum of its Fourier 
coefficients within the bounds of Equation 3.3 such that 


Pwin{{ x i V)'i t) 


00 oo 

12 12 Anm ex p 

771= —oo n= — oo 


m mx\ 
— ) exp 


inny 

Ly 
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where A nm are the corresponding Fourier coefficients. 

Since p((x, y);t) must be symmetric due to the nature of its construction, then A- n ~ m = 
A nm • The individual terms of Equation 3.8 can now be written as 
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Expanding this and rewriting the whole summation from Equation 3.8 gives 
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(3.10) 


Because the wire locations on the wave gage are fixed, the vector between any two gage 
elements has a known value ( x p ,y p ) for p = 1,2,3,... ,p m ax ■ The covariance function 
evaluated at the location of each of these vectors is approximated as the time average of 
the products of the surface elevations at the two corresponding gage elements, as per its 
definition. Since only a finite set of measured surface elevation data is available, discrete 
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values of the covariance function, p((x p ,y p )-,t), can be calculated from the finite set of 
measured surface elevation data. The series from Equation 3.10 can then be approximated 
by its truncated form: 


Pwin{{%pi Up) i — ^00 + 2 'y ] A n o cos 


+ 2£ £ + (3.11) 

m=ln=-(N-l) V L * L V ) 


The unknowns in this equation are the Fourier coefficients, which are the A nm s. There are 
M(2N — 1) — IV + 1 of these unknown coefficients. If enough (well behaved) values of the 
covariance are known, this system of equations can be solved for the A nm s. In other words, 
if 


Pmax > M(2N - 1) - TV + 1 (3.12) 

then one can use a minimum mean square error matrix solver technique to find the Fourier 
coefficients of the covariance matrix by solving a matrix equation of the form 

Pwin{x p ,y p j\ = [Aim] c COs(^^ £ + -y^) (3.13) 

j Lj X Jjy 

where C is a set of constant coefficients determined by putting Equation 3.11 in this matrix 
form. 


The Directional Spectrum 

The classical method of determining a two-dimensional spectrum of a spatially two-dimensional 
signal is to take the two-dimensional Fourier transform of the covariance of that signal. For 
a continuous signal this corresponds to computing the following integral 

[ r oo roc 

$(k x , ky) = ?—-T 2 / / p(x, y) exp [-i (k x x + k y y)] dxdy (3.14) 

{Z7T) J _oo J —oo 
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If p(x , y) is replaced by its representation from Equation 3.8 and then the exponential terms 
are simplified, the following form for the two-dimensional spectrum results 


OO CO 
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The integrals can now be separated and individually evaluated giving 
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Simplifying this slightly gives the final result of 

OO OO 
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(3.17) 


This equation can be used to calculate the spectrum at any wave number values, (k x ,k y ), 
assuming the A nm coefficients are known. For very particular values of k x and k y the 
formula above simplifies even further. This is for: 
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From Equation 3.11 the A nm coefficients can be calculated for m = 1,2,... , M — 1 and 
n = -{N - 1), -(IV - 1) + 1,... ,0,... ,(N- 1). 

The equations above provide the key analytical tools for both estimating the covari¬ 
ance function given a finite set of surface height measurements, and evaluating the two- 
dimensional spectrum at either specific wavenumber locations or at any location, given a 
Fourier coefficient expansion of the covariance function. The limitations of this method 
stem primarily from the accuracy of the covariance estimate in Equation 3.11; if there are 
not enough values of the covariance function or the minimum mean square error fit is not 
well behaved, the resulting Fourier coefficients will not be good estimates, leading to inac¬ 
curate results using Equations 3.17 through 3.20. Every effort has been made to assure that 
these errors are as small as possible. Some tests to determine the accuracy of the analysis 
methods were described in the author’s Masters thesis [61]. 

Calibration 

If both the input and output of a linear system are known in time, they can be transformed 
into the frequency domain and used to generate the frequency response function as described 
in Section 3.3.2. Once the frequency response has been calculated, it can be used with the 
system output to reconstruct the input signal - this is how the signal outputs from the wave 
gage can be used to regenerate the measured wave heights. The goal of calibration is to 
produce the frequency response. 

For any analog electronic system like the wave gage, the shape of the frequency response 
function can be modeled using the nominal values of the various circuit elements that 
make up the frequency-dependent portions of the analog filtering. This was illustrated 
in Table 3.1. This model can predict the shape of the frequency response, but it cannot 
accurately predict the magnitude or the exact location of the peak of the frequency response. 
The reason for this is that the values of the circuit elements have some uncertainty associated 
with them. For instance, capacitors are generally purchased with specific tolerances such 
as 5 or 10%. The exact values of all the circuit elements are not known unless measured 
independently. More important, the values of these elements can and do change with 
temperature. The modeled theoretical frequency response cannot, therefore, be used alone 
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to represent the precise characteristics of the system. 

Tests were conducted to measure the temperature dependency by putting the analog 
electronics in a temperature-controlled oven. The resulting affects on the calibration of each 
wire as a function of temperature are used during the calibration of the instrument if there 
is a significant temperature change over the course of the experiment. 

The frequency response of the filtering system can most easily be measured by ta kin g 
the Fourier transform of the time series output of the system’s response to an impulse input 
signal. Unfortunately, generating a ‘real-world’ impulse response is rather difficult, and due 
to the nature of the circuitry, it would be difficult to insert the impulse into the filter input 
if one could be generated. Generating a step input in the system by stepping the amplitude 
of the high-frequency carrier wave is much easier to accomplish and contains the relevant 
information that an impulse response would contain. This relationship between Fourier 
transform of the step response, F step (w), and the transfer function, is 

H(u>) = Y step (ui) • ^ (3.21) 

where h is the effective length in centimeters of nichrome wire that is shorted out to generate 
the step which is a predetermined constant. 

Now that the frequency response for each channel has been both measured and modeled, 
all that remains is to fit the model to the actual measurement to establish a ‘best-fit’ model. 
This new model will be used to ‘unfilter’ the output signal from each of the channels 
and return the wave height. The reason a fitted model must be used instead of simply 
using the actual measured frequency response is because the measured frequency response 
only contains magnitude information, but the actual frequency response contains phase 
information as well. The phase information is tied up in the temporal relationship between 
the step and its response. This temporal coupling is not well measured. A fitting routine 
is used instead which adjusts the nominal values for the capacitors and resistors in the 
transfer function model until a best fit between the resulting shape and the measured shape 
is found. 
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3.7.2 Program Implementation 

Each analysis step in the data processing package section of Figure 3-1 is performed by a 
separate program. This was done because Matlab, FORTRAN, and C coding are each used 
to accomplish various analytical operations, to facilitate debugging and testing, and because 
the programs are each very lengthy both in file size and in processing time. The following 
sections contain brief descriptions of each of the analysis programs. These descriptions are 
not meant to expand on the discussion of the signal processing and calibration methods 
presented earlier in this Chapter and in Chapter 2, but are instead meant as a guide both 
to the use of the code and to the implementation of the theory. 

Calibration (step 1): calib.m 

Input Files: stepl.asc, step2.asc, step3.asc, step4.asc, 

zerol.asc, zero2.asc, zero3.asc, zero4.asc, and <run 
numeric>. asc 

Prompted Input: <run numeric> (the file name of the data file to be calibrated) 
Output: data.out and badchan.out (input to postcal.for) 

To calibrate the wave gage buoy output a fairly detailed set of steps must be followed: 

1. Convert the <run numeric> .wav file into an ASCII file using wav2asci. c. 

2. Convert the 8 calibration files from binary into ASCII files using wav2asci. c. 

3. Rename the files containing step information to step 1. asc, step2.asc, step3.asc, 
and step4.asc, and rename the files containing noise information to zerol.asc, 
zero2.asc, zero3.asc, and zero4.asc. 

4. Copy the step* and zero* files into a directory with the calib.m file. 

5. Open Matlab in this directory and run the calib.m Matlab file. 
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3.7.3 Calibration (step 2): postcal.for 

Input Files: data-out, badchan.out (output from calib.m), and 

<run numeric>. KVH 

Prompted Input: <run numeric> (the file name of the data file to be calibrated) 

Output: direct, out, calstats. out, stats, out, 

ch*spec.out, mssstats.dat, newchan*.out, and autoin.dat 
(input to sw_auto. c) 

These programs calibrate the wave wire voltage output and return wave height output 
at each wire location as a function of time. This is done by comparing the measured 
impulse response of each of the channels with a theoretical one generated as was described 
in Section 3.7. The theoretical impulse response is then ‘tuned’ so that it best matches the 
calibration data collected. The uncalibrated wave wire data is then divided by the ‘tuned’ 
impulse response in the frequency domain for every channel and then it is converted back 
into the time domain. The resulting time series have effectively had the analog pre-whitening 
filter backed out of them. 

This calibration program also incorporates the compass information by reassigning the 
channel 1 time series to be the time series of the wire which is closest to magnetic North at 
each scan. Each of the other channels is reassigned in a corresponding manner relative to 
channel 1. This allows an absolute directional spectrum to be calculated relative to magnetic 
North. There is a little over 5 degrees of error in the resulting directional calculations due 
to the 5.6 degree spacing between each of the wires. 

Finally, if any of the wires is determined to be broken or behaving badly, the program 
interpolates the wave heights between the two adjacent channels. If more than two adjacent 
channels are determined to be broken the data set is not used. This was not a problem for 
the data sets presented later in this thesis; only 4 channels needed to be interpolated. 

3.7.4 Calculation of Covariance: correl.c 

Input Files: autoin.dat (output of autocal*.for) 

Prompted Input: number of scans to evaluate 

Output: correl.mas (input to dfsum*.for) 
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This program calculates the covariance of the sea-surface elevation data. This is done 
by multiplying each data point of a particular scan by every other data point in that scan. 
Each of these covariance values is a function of the difference in locations of each of the 
data points used to generate it. The ensemble average of the covariance of each of the 4096 
scans is then calculated and printed to a file. Unlike a covariance calculated over time or 
in one dimension, the covariance calculated in this program is a function of translation in 
two dimensions, x and y. A more explicit description of this program is contained in Chen’s 
Masters thesis[12] where it is called sw_auto.c. 

3.7.5 Calculation of Directional Spectrum: dfsum*.for 

Input files: correl.mas (output of correl.c) 

Prompted Input: none 

Output Files: s.tec, z.tec, sk.tec, anglspec.out, ptspec.out 

This program calculates the directional spectrum of an input data set which contains 
the magnitude of the covariance of the data set and the location at which these covariance 
values are known. First the covariance data set is convolved with a smoothing window, 
then the truncated Fourier series coefficients are calculated for the covariance data set, and 
finally, these Fourier coefficients are used to estimate the directional spectrum at particu¬ 
lar wavenumbers in both the North-South and East-West direction. This method is fully 
described in the author’s Masters thesis [61]. 

In addition to generating the directional spectrum, this program also produces an an¬ 
gle spectrum and a point spectrum. The angle spectrum is a representation of the spec¬ 
trum along various slices taken along a particular radial direction from the origin. The 
point spectrum is calculated by integrating the two-dimensional spectrum over a particular 
wavenumber. These are used to compare the results of the spectral calculations to other 
experimentally or theoretically generated spectra data collected by other researchers. 
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Chapter 4 


Experimental Results 


The at-sea experiment which produced the bulk of the data presented in this document was 
conducted on October 28, 1998 with the aid of Sam Tolkoff, Monika Bang-Campbell, and 
Dave Olmstead, the captain of the Woods Hole Oceanographic Institution’s R/V Asterias. 

4.1 Buzzards Bay Tower Experiment 

The Buzzards Bay Tower (BUZM3) is located in the southwest entrance to Buzzards Bay, 
east of Newport, Rhode Island at 41.40N latitude, 71.03W longitude. It is a triangular 
C-MAN station which is monitored by the National Oceanographic and Atmospheric Ad¬ 
ministration (NOAA) and the Woods Hole Coast Guard Aids to Navigation Team. The 
tower is approximately 25 meters high with two platforms — the lower platform is approxi¬ 
mately 2 meters above the mean high water mark with a lower docking level to allow easy 
access by boat, and the upper platform is approximately 15 meters above mean sea level 
and is where the bulk of the NOAA instrumentation is located as illustrated in Figure 4-1. 

4.1.1 Experimental Setup 

The instruments were mounted on the lower platform in the positions depicted in Figures 4-2 
and 4-3. The circular wire wave gage buoy was deployed just East of the tower. 


62 




63 











Figure 4-2: BUZM3 Instrument Layout, Top View (not to scale) 



Figure 4-3: BUZM3 Instrument Layout, Side View (not to scale) 
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Circular Wire Wave Gage 


The circular wire wave gage was designed to float freely with the surface waves and currents 
to reduce flow distortion at the wire locations central to the buoy frame. For this experiment 
it was deployed outside of the expected wind and wave shadow of the tower and allowed 
to drift for the duration of the experiment. The R/V Asterias occasionally approached the 
buoy from down wind to obtain a Global Positioning System (GPS) fix. The drift path of 
the buoy is shown in Figure 4-4. Because the wind was blowing from the open ocean out 
of the southwest it is assumed that the conditions at the buoy were statistically similar to 
those experienced at the tower and that the fetch was over 100 kilometers. The similarity 
of surface wave conditions is important for the comparisons made between measurements 
collected at the tower and those collected by the buoy. The buoy drifted approximately 
2.5 nautical miles over the course of its deployment, however, and the validity of even 
meteorological similar conditions between the buoy and the tower at these distances could 
be cited as a possible reason for the discrepancies between the resulting measurements at 
the two locations. 

The buoy was deployed at 1220 EST and was preprogrammed to collect 6 data sets of 
293 blocks each (~ 25 minutes) at 30 minute intervals beginning at 1215 EST. The buoy 
was recovered at 1506 EST. The first and last data sets were not used because the buoy 
was not in the water and free of the deployment vessel for their entirety. 

Scatterometer 

The scatterometer is a dual-polarized continuous-wave X-band Doppler backscatterometer 
designed and built at the Johns Hopkins University Applied Physics Laboratory and was 
on loan to WHOI. It has four antennas, two of which transmit the H- and V-polarized fields 
at 10.00 GHz and 10.06 GHz, respectively, and two which receive the backscattered signals. 
The output from the scatterometer is a time series of in-phase and quadrature signals in 
both HH and VV polarizations (i.e. an I and a Q output for each polarization). The output 
signals were digitized at 1000 Hz. 

The scatterometer was mounted on an L-shaped frame which was bolted to a stand 
just above the railing of the lower platform of the tower. The frame extended horizontally 
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Figure 4-4: Circular Wire Wave Gage Buoy Drift Track. The Buzzards Bay Tower is located 
in the bottom left corner of the map. The gray dots with the thick line connecting them 
beginning just to the East of the Tower represent the drift track of the buoy. The gray 
dots represent the location fixes obtained by the Asterias during deployment. For reference 
purposes True North is vertically upwards and the entire image is approximately 5.5 nautical 
miles wide in the E-W direction. The depths given on the chart are in feet. 
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outward 1.47 meters from the railing and 1.27 meters above the railing. The railing was 
approximately 1 meter high. Three fiber lines were used to further secure the frame - two 
extended to either side of the frame and were tied to the railing, the third was secured 
vertically above the frame. The frame also allowed the scatterometer to be mounted at a 
user specified angle of declination - an incidence angle of 59° was used for this experiment. 
Additionally, the very top portion of the frame housed a motor-driven rotary mechanism 
which pivoted the scatterometer in an arc of more than 180°, allowing a field of view 
encompassing the entire windward face of the tower. The scatterometer completed one 
180°+ rotation in approximately 28 seconds. At each end of the rotation a limit switch 
was triggered which reversed the direction of the rotation. An additional 15° of rotation at 
each end of the arc allowed enough extra time for the small motion caused by the abrupt 
change in direction of the scatterometer to be adequately damped, so as to not influence the 
scatterometer’s measurements during the 180° arc when it was actually viewing the open 
water and not the tower. 

At 1346 EST on the day of the experiment, the mean water level at the tower was 
approximately 2 meters below the grating on the lower platform. With this water height and 
the dimensions of the frame given above, the center of the field of view of the scatterometer 
was just over 7 meters from the tower and the approximate Bragg wavelength was 1.7 
centimeters. The distance of the scatterometer footprint was maximized to ensure that any 
waves reflected from the tower influenced the scattered return as little as possible. 

The scatterometer data shown later in this chapter is not presented in calibrated ctq 
form. Because the raw backscattered power is directly proportional to ctq through a constant 
calibration factor, and the only aspect of the scattered return of significance to this study is 
its azimuthal shape (a purely relative quantity), it was not necessary to absolutely calibrate 
the instrument. 

Capacitive Wire Wave Gage 

A capacitive wave gage array with six wires was hung from the frame on which the scat¬ 
terometer was mounted. The capacitive wires were constructed following the guidelines 
described by Chapman and Monaldo [10], the electronics were designed by Andrey Karach- 
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intsev at the University of Rhode Island (URI), and the processing code used to generate 
two dimensional spectra was written by Kurt Hanson [23], also from URI. The instrument 
operates by measuring the capacitance between the insulated core of the wave wire, and 
the conductive salt water surrounding it. Each instrument outputs voltages that are digi¬ 
tized at 1000 Hz. Calibration measurements are then used to relate the voltage level to the 
appropriate capacitance value. The calibration tables for the six wires used are given in 
Appendix I. The resulting capacitance values can be related to the level of the salt water 
surrounding the wire and, thus, be used to determine sea surface height. 

The wires were distributed unevenly around a hexagonal frame with a minimum spacing 
between wires of 0.3 meters and a maximum spacing of 0.97 meters as pictured in Figure 4- 
5. The upper and lower hexagonal frames which supported the electronics and 1.6 meter 
long wires were separated by a 1.8 meter long thin vertical pole located in the center of 
the hexagon. Changes in the mean water level due to the tides required that the frame 
occasionally be moved vertically in order to insure that the measured waves not be clipped 
by limited vertical range of the wires. 
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Figure 4-5: Capacitive Wire Wave Gage Hexagonal Frame and Wire Location Diagram. 
Dimensions are in inches. The hexagonal frame is shown with a dashed line and the wire 
locations are indicated with open circles 
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The top of the array was secured to the scatterometer frame with a pulley and rope and 
was also tied to the railings to minimize twisting. Because the lower hexagonal frame was 
not attached to anything, it tilted and rotated slightly with the passing of the larger waves. 
It was estimated by observation that errors caused by these unavoidable motions were not 
greater than 2 centimeters in wave height and 10° in angle. 

Sonic Anemometer 

A CSAT3 Campbell Scientific Inc. (Logan, Utah) three-dimensional sonic anemometer was 
mounted on a triangular antenna frame and extended about 2.6 meters away from the 
tower base. The boom was extended parallel to the face on which the scatterometer was 
mounted, and the sensor was turned into the wind to allow for minimum flow distortion as 
illustrated in Figure 4-2. The location of the sensor was approximately 2.6 meters above 
the mean water surface. The boom was secured from above with a line which, along with 
the bolting of the boom to the tower platform, minimized the motion of the sensor. The 
analog outputs from the sonic anemometer were calibrated using the low range calibration 
constants provided by the manufacturer in the Campbell Scientific CSAT3 Manual, page 7, 
and shown in Table 4.1. 

Ux Uy Vjz 

Slope (m s" 1 V- 1 ) : 6.5536 6.5536 1.6384 
Table 4.1: Sonic Anemometer Calibration Coefficients 


4.1.2 Buzzards Bay Tower NOAA Data 

The NOAA instrumentation on the BUZM3 C-MAN tower includes the following measure¬ 
ments: wind direction, wind speed, wind gust magnitude, air temperature, and barometric 
pressure. The anemometer is located 24.8 meters above mean sea level, the air temperature 
sensor is located 24.5 meters above mean sea level, and the barometer is located 17.4 meters 
above mean sea level. The other sensors at the tower were either broken or not reliable as 
recorded at the National Data Buoy Center station information web page on October 28th, 


69 





1998. For further information on the tower and its instrument suite, please reference either 
the National Weather Service web site http://www.nws.fsu.edu or the National Data Buoy 
Center web site http://www.ndbc.noaa.gov. Both web sites contain similar information and 
allow access to archived data. 

4.1.3 Results 

General Ocean and Atmospheric Conditions 

The atmospheric conditions during the experiment on October 28th were fluctuating. The 
barometer was dropping, the air temperature was rising, and the wind velocity was slowly 
increasing. These trends can be seen in Figures 4-6 and 4-7. 1 They are shown for the sole 
purpose of establishing the synoptic conditions during the experiment. The more detailed 
wind velocity and direction information shown in Figures 4-8, 4-9, and 4-10 2 is used for the 
spectral model input parameters used in Chapter 5. A histogram showing the distribution 
of the wind direction measurements (unaveraged) during the first data set (10281315) is 
also included in Figure 4-11. Histograms of the other two wind direction data sets show 
a similar shape. Over the course of the experiment the location of the peak of the wave 
height frequency spectrum remained constant at approximately 0.30 Hz as measured by 
time series spectral calculations using the capacitive wire wave gage data. 

The long wave properties were measured using the capacitive wire wave gage. Each wire 
was sampled at 1000 Hz, but the instrument’s maximum resolution is only about 4 Hz due 
to electronic and wire limitations, so the data was subsampled down to 100 Hz. This lower 
frequency is still oversampled, but it prevented any high frequency signal or noise aliasing 
which can be verified by comparing a single channel frequency spectrum from both the 
original and the subsampled data sets shown in Figure 4-12. Figure 4-13 shows an average 
of the one-dimensional frequency spectra from each of the six wires for each data set. These 
three plots show a consistency both in the location of the peak of the spectrum as well as 

1 All azimuthal experimental measurements in this thesis are given in Magnetic North coordinates, except 
where noted, and will be referred to simply as North. At the time of this experiment the variation between 
Magnetic North and True North was 15° 31’ W. 

2 The wind velocity data shown in Figures 4-8, 4-9, and 4-10 was adjusted from its measurement height of 
2.6 meters to a standard 10 meter height using the Charnock model and a neutral stability assumption [11], 
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Figure 4-6: Hourly Barometric Pressure and Air Temperature from BUZM3 Tower. The 
barometer is located 17.4 meters above mean sea level and the temperature sensor is located 
24.5 meters above mean sea level. The data points are marked with stars. The dashed lines 
indicate the approximate temporal duration of the experimental data collected with the 
author’s instruments. 


71 










Local Time (24 hour clock) 



Local Time (24 hour clock) 

Figure 4-7: 10 Minute Averaged Wind Direction and Velocity from BUZM3 Anemometer. 
The anemometer is located 24.8 meters above mean sea level. The data points are marked 
with an ‘x’. The dashed lines indicate the approximate temporal duration of the experi¬ 
mental data collected with the author’s instruments. The squares indicate the maximum 
wind velocity during each hour and when it occurred. 
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13.3 13.4 13.5 13.6 13.7 



13.3 13.4 13.5 13.6 13.7 

Time (24 hour clock) 

Figure 4-8: CSAT3 Sonic Anemometer Data from Data Set 10281315. The wind velocity is 
given at a 10 meter reference height. The data has been averaged over 2 second intervals. 
The mean wind direction is 213.7° and the mean wind velocity is 8.05 m/s. 
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Figure 4 9. CSAT3 Sonic Anemometer Data from Data Set 10281344. The wind velocity is 
given at a 10 meter reference height. The data has been averaged over 2 second intervals. 
The mean wind direction is 220.9° and the mean wind velocity is 8.23 m/s. 
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Direction (degrees N) 



Time (24 hour clock) 

Figure 4-10: CSAT3 Sonic Anemometer Data from Data Set 10281414. The wind velocity 
is given at a 10 meter reference height. The data has been averaged over 2 second intervals. 
The mean wind direction is 225.6° and the mean wind velocity is 7.96 m/s. 
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Figure 4-11: Wind Direction Histogram for Data Set 10281315. 



Figure 4-12: Capacitive Wire Wave Gage Single Wire Frequency Spectrum from Data Set 
10281344. The pluses are the spectrum calculated from the first 8.7 minutes of the data set 
which was sampled at 1000 Hz. The solid line represents a spectrum calculated from the first 
10.9 minutes of the data set which was a subsampled version of the original with a sampling 
rate of 100 Hz. The spectrum is shown only out to the resolution of the instrument. The 
noise floor for each channel plateaus at around 10 Hz between 10“ 2 and 10 -3 cm 2 /Hz. 
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(a) (b) (c) 


Figure 4-13: Capacitive Wire Wave Gage One-Dimensional Frequency Spectra Averaged 
From All 6 Wires From Data Sets 10281315, 10281344, and 10281414. Subfigure (a) was 
calculated from the first 21.8 minutes of data set 10281315, while subfigures (b) and (c) were 
calculated from the first 10.9 minutes of data sets 10281344 and 10281414, respectively. 


its shape indicating that while the atmospheric conditions were not stationary, the change 
in the long wave field energy density between the first and last data sets was not significant. 

The directionality of the long wave spectrum is also available from the capacitive wire 
wave gage array. Plots of the two-dimensional frequency spectrum from each of the data 
files are given in Figures 4-14, 4-15, and 4-16. The direction from which the wind is blowing 
is indicated on each plot and correspond well with the calculated energy distribution which 
should appear ‘downwind’ on these plots due to the processing algorithm used (e.g. in Fig¬ 
ure 4-14 the wind is coming from approximately 214° and the energy is bunched near the 
‘downwind’ direction of 214-180=34°). Their purpose, however, is to show general direc¬ 
tional energy distribution, not specific details, as the resolution and method of deployment 
for the instrument were not designed for this. Additionally, the low frequency limit in the 
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figures of 0.35 Hz is the result of the directional resolution limit warranted by the maximum 
wire spacing in the array: a 0.35 Hz wave has a wavelength of approximately 12.6 meters. 
The high frequency directional limit of the instrument is 2.29 Hz corresponding to a wave 
having a 0.3 meter wavelength. The plots do not extend out to this higher limit because 
the energy density is plotted on a linear scale, and the power law fall-off of energy density 
would prevent the representation of energy levels at these higher frequencies. The peaks 
in the directional energy distribution plots between 0.4 and 0.5 Hz are secondary peaks in 
the spectrum which can be seen in the higher resolution point spectrum subfigure (a) in 
Figure 4-13. 

These spectra are computed using a two-dimensional form of the Data Adaptive Spectral 
Estimator (DASE) proposed by Davis and Regier [15]. This method attempts to estimate 
a continuous two-dimensional spectrum, $(&), where k is a vector, using the following 
approximation 


< $(fc) >= J dk'$(k')W(k,k')W(k,k') = £> nm (fc)exp(ifc' ■ [x n - x m )) 

n,m 

(4.1) 


where W(k,k') is the 


wavenumber window and measures how much an estimate of the spectrum 
at wavenumber k is affected by the variance density at wavenumber k' in ex¬ 
actly the way that the sine function frequency window encountered in conven¬ 
tional time series analysis measures variance spreading in frequency space. Data 
adaptive spectral estimators do not try to approximate W with a delta func¬ 
tion, however, as in conventional a priori analysis methods. Instead the DASE 
method attempts to make W(k, k') small at those values of k' ^ k where <f>(fc) 
is significant. 

- Davis and Regier [15] 

Another noteworthy detail about this analysis method is that it constrains the wave energy 
in the estimated spectrum to lie only along the dispersion curve. While this limitation is 
an approximation, as discussed in Section 2.1, for gravity waves it is a reasonable one and 
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Figure 4-14: Capacitive Wire Wave Gage Two-Dimensional Frequency Spectrum from Data 
Set 10281315. The spectral estimate has units of cm 2 /(Hz 2 ) and gets darker as the energy 
level increases. North is pointing toward the top of the page and East is pointing to the 
right. The dashed line from the edge into the center of the plot shows the mean wind 
direction during the data set. Energy is only shown between the bands of 0.35 Hz and 1 
Hz. This estimate is based on the first 21.8 minutes of the data set. 
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Figure 4-15: Capacitive Wire Wave Gage Two-Dimensional Frequency Spectrum from Data 
Set 10281344. The spectral estimate has units of cm 2 /(Hz 2 ) and gets darker as the energy 
level increases. North is pointing toward the top of the page and East is pointing to the 
right. The dashed line from the edge into the center of the plot shows the mean wind 
direction during the data set. Energy is only shown between the bands of 0.35 Hz and 1 
Hz. This estimate is based on the first 10.9 minutes of the data set. 
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Figure 4-16: Capacitive Wire Wave Gage Two-Dimensional Frequency Spectrum from Data 
Set 10281414. The spectral estimate has units of cm 2 /(Hz 2 ) and gets darker as the energy 
level increases. North is pointing toward the top of the page and East is pointing to the 
right. The dashed line from the edge into the center of the plot shows the mean wind 
direction during the data set. Energy is only shown between the bands of 0.35 Hz and 1 
Hz. This estimate is based on the first 10.9 minutes of the data set. 
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is widely used for long-wave spectral estimation formulations. 


Scatterometer Measurements 

The scatterometer data was processed as follows. The raw output files constitute a time 
series of voltages collected as the instrument rotates 180°+ clockwise and then 180°+ coun¬ 
terclockwise. First the raw data was separated into periods corresponding to specific sectors 
of rotation. The angular rotation was measured using a potentiometer mounted on the mo¬ 
tor axis. The power for each polarization was calculated for each sector by squaring the 
in-phase and quadrature measured signals and then adding them together. Each consecu¬ 
tive arc traveled the opposite direction from the previous one, so the power from alternate 
data periods was averaged. Then each was shifted so the same geographic directional sec¬ 
tors were aligned and then they were averaged together. The resulting directional power 
averages, one for HH polarization and one for VV, were referenced to magnetic North by 
comparing the known rotation rate and total rotation angle of the scatterometer with its 
position when mounted on the tower. The uncertainty in absolute direction was about 5°. 
It should be noted that the time between each illumination by the scatterometer of a specific 
location on the water surface was not equal. This was due to the reversal of direction of the 
scatterometer scan. Even though the ocean surface is not truly ergodic either on the day 
of the experiment or in general, the uneven spacing of the data samples which are averaged 
is assumed not to affect the resulting directional power plots. 

Figures 4-17, 4-18, and 4-19 show the averaged backscattered power measured by the 
scatterometer for both polarizations during each of the data sets. The total averaging 
time for each data set is approximately 26 minutes. The data represents only the 180° of 
rotation when the scatterometer was viewing open water (the extra 15° at each end of the 
rotation has been removed for clarity). The data was not entirely circular, meaning there 
may have been small variations in the phase shift of the output signals, so corrections were 
made including a general correction to equilibrate overall gain between the I and Q channels 
for each polarization. The description of these modifications can be found in an internal 
report by J. R. Jensen [35] at the Johns Hopkins University Applied Physics Laboratory. 
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Figure 4-17: Directional X-band Scatterometer Return from Data Set 10281315. The power 
is calculated using I 2 +Q 2 where I and Q are measured in volts. The thick solid line rep¬ 
resents a best fit of the data to the function cos 2s while the thick dashed line 

represents a best fit of the data to the function cosh~ 2 (/?(0 — 0 max )). The s and (3 values 
for the HH and W polarizations are 4.02, 0.879 and 5.31, 0.890, respectively. The data 
points shown are 85 point averages. The thin vertical dashed line is the mean wind direction 
during the data set of (214°) 
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Figure 4-18: Directional X-band Scatterometer Return from Data Set 10281344. The power 
is calculated using I 2 +Q 2 where I and Q are measured in volts. The thick solid line rep¬ 
resents a best fit of the data to the function cos 2s while the thick dashed line 

represents a best fit of the data to the function cosh -2 (/?(<? — 9 max )) • The s and (3 values 
for the HH and W polarizations are 3.79, 0.883 and 5.44, 0.878, respectively. The data 
points shown are 85 point averages. The thin vertical dashed line is the mean wind direction 
during the data set of (221°). 
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Figure 4-19: Directional X-band Scatterometer Return from Data Set 10281414. The power 
is calculated using I 2 +Q 2 where I and Q are measured in volts. The thick solid line rep¬ 
resents a best fit of the data to the function cos 2,s j, while the thick dashed line 

represents a best fit of the data to the function cosh ~ 2 {f3(8 — 0 rnax )). The s and f3 values 
for the HH and VV polarizations are 3.72, 0.831 and 5.73, 0.839, respectively. The data 
points shown are 53 point averages. The thin vertical dashed line is the mean wind direction 
during the data set of (226°). 
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Figure 4-20: Circular Wire Wave Gage Frequency Spectrum Averaged from 16 Single Wire 
Spectra from Data Sets 10281415, 10281445, and 10281515. Subfigure (a) is from data set 
10281415, (b) is from data set 10281445, and (c) is from data set 10281515. 


Circular Wire Wave Gage Measurements 

The frequency spectra from approximately 21 minutes of calibrated time series data from the 
three different circular wire wave gage data sets are shown in Figure 4-20. Each spectrum 
is an average of the frequency spectra from the first sixteen channels of the instrument. 
All the data set file names are based on the month, day, hour, and minute of collection 
corresponding to Eastern Daylight Time (EDT), however, the circular wire wave gage used 
Eastern Standard Time (EST, where EST-1 = EDT), so while the data sets were collected 
at 1315, 1345, and 1415 EDT, their file names are 10281415, 10281445, and 10281515, 
respectively. 

Two-dimensional wave height spectra are presented in Figures 4-21, 4-22, and 4-23, 
where each spectrum was calculated from data averaged over approximately 21 minutes. 
It should be noted that the compass directions on these plots are not standard: North is 
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upwards on the page, but West is to the right. Additionally, the spectrum is symmetric 
through the origin because of its definition as being the Fourier transform of an autocorre¬ 
lation which is, itself, symmetric through the origin. Because of this symmetry only half of 
the spectrum is shown in the two-dimensional spectrum plots. 

A one-dimensional wavenumber spectrum can be generated from each of the two-dimensional 
spectral plots by integrating out the angular dimensions, such that given a symmetric two- 
dimensional spectrum of the form <h(/c, 9) with 0 only given over a 180° range, (Figures 4-21 
through 4-23 can be converted to this form), then an omnidirectional spectrum can be 
generated using 

S(k) = 2 r ${k,9)kd9 (4.2) 

Jo 

The omnidirectional spectra calculated in this manner from the circular wire wave gage 
two-dimensional wavenumber spectra are plotted in Figure 4-24. 

A more detailed look at the directional spectrum plots can be obtained by taking circular 
slices of the spectrum. Figure 4-25 shows some of these slices, each depicting the actual 
directional distribution of the energy at a given wavelength. The white areas in the two- 
dimensional wavenumber spectral plots and the points along the directional energy slices 
which are equal to zero represent areas where the energy estimate is either zero or negative. 
The processing algorithms allow negative energy due to small negative side lobes caused by 
filtering, but negative energy is not physically consistent, so any negative energy areas were 
set to zero. In general, the magnitude of the negative energy lies within the magnitude of 
the uncertainty of the spectral estimate. 

4.2 Comparison of Data Processing Methods 

In order to verify that the analysis method being used by the circular wire wave gage does 
indeed reproduce the appropriate estimate of the spectrum, it was compared to the data 
processing used by a scanning laser slope gage (SLSG) which also has the ability to measure 
the two-dimensional wavenumber frequency. If an SLSG calibrated slope time series data 
is used to estimate the slope spectrum using the circular wire wave gage processing, and 
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Figure 4 21: Circular Wire Wave Gage Two-dimensional Wave Height Spectrum from Data 
Set 10281415 (corresponding to 1315 on October 28, 1998). The spectral energy, $(fc x , k y ), 
has units of cm 2 /(rad/cm) 2 and the colorbar scale reflects the spectral energy plotted as 
10 (> ky))- The darker the color, the higher the energy level. North is pointing toward 
the top of the page and West is pointing to the right. The dashed line is the mean wind 
direction as measured at the BUZM3 tower over the approximate time period that this data 
set was collected. This estimate is based on the first 21 minutes of the data set. 
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Figure 4-22: Circular Wire Wave Gage Two-dimensional Wave Height Spectrum from Data 
Set 10281445 (corresponding to 1345 on October 28, 1998). The spectral energy, Q(k x ,k y ), 
has units of cm 2 /(rad/cm) 2 and the colorbar scale reflects the spectral energy plotted as 
log 10 (4>(fc x , k y )). The darker the color, the higher the energy level. North is pointing toward 
the top of the page and West is pointing to the right. The dashed line is the mean wind 
direction as measured at the BUZM3 tower over the approximate time period that this data 
set was collected. This estimate is based on the first 21 minutes of the data set. 
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Figure 4 23: Circular Wire Wave Gage Two-dimensional Wave Height Spectrum from Data 
Set 10281515 (corresponding to 1415 on October 28, 1998). The spectral energy, $(fc x , k y ), 
has units of crn^/(rad/cm) ^ and the colorbar scale reflects the spectral energy plotted as 
ky))- The darker the color, the higher the energy level. North is pointing toward 
the top of the page and West is pointing to the right. The dashed line is the mean wind 
direction as measured at the BUZM3 tower over the approximate time period that this data 
set was collected. This estimate is based on the first 21 minutes of the data set. 
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Figure 4-24: Omnidirectional Spectra from the Circular Wire Wave Gage Data Sets 
10281415, 10281445, and 10281515. Each spectrum was generated using 21 minutes of 
data. Subfigure (a) is from data set 10281415, (b) is from data set 10281445, and (c) is 
from data set 10281515. 
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Figure 4-25. Directional Slices of Circular Wire Wave Gage Two-Dimensional Wavenumber 
Spectra from Data Sets 10281415, 10281445, and 10281515. The axes on subfigures (b) - 
(d) are the same units as those in (a). Note that the y-axis scale of subfigure (b) is not the 
same as the others, however, and the x-axis has been shifted for (b) and (d). The numbers 
next to each curve represent the magnitude of the wavenumber at which the slice was taken 
and are given in rad/cm. Subfigures (a) and (b) are from data set 10281415, (c) is from 
data set 10281445, and (d) is from data set 10281515. 


92 








the shape of this spectrum compares well with the shape of the slope spectrum estimated 
using the SLSG processing, then it can be said that the two methods, given the same input 
data, generate the same or similar estimates of the spectrum, thus validating the CWWG 
processing methods. 

The Scanning Laser Slope Gage was developed in the early 90’s by Dr. Erik Bock at 
the Woods Hole Oceanographic Institution and Dr. Tetsu Hara at the University of Rhode 
Island to measure the three-dimensional wavenumber-frequency slope spectrum of capillary- 
gravity waves. It was designed to be deployed from a catamaran-type platform with two 
subassemblies: a laser pod which produces a collimated laser beam which is rotated in a 
circle producing a circular scan, and a head unit that measures the angular deflection of 
the incoming light beam using an array of photodiodes. The laser pod is located under the 
water with the beam centerline facing the air-water interface, and the head unit is lined 
up vertically above it. Both are mounted forward of the bow of the deployment platform. 
When the beam from the laser pod breaches the air-water interface it is refracted according 
to Snell’s law. The angular deflection of the refracted beam is used to determine the surface 
slope at the point the beam left the water surface relative to the position of the SLSG. The 
result is a slope measurement at each successive position of the laser beam around the circle 
it traces. For the data presented in this section each scan was collected at 138.8 Hz. 

The three-dimensional spectrum is estimated by calculating the Fourier transform of 
the autocorrelation function, or covariance, of the surface slope. Using the same definitions 
as were used in equation 2.3, the sea surface slope can be represented by 

dy dy 

— = si{x,y,t) — = s 2 {x,y,t ) (4.3) 

where y = (sc, t ) = (sc, y, t ) is the sea surface displacement. The autocorrelation function of 
the zth component of the surface slope can then be defined as 


R?(Z,M = 8i{x,y,t)si(x + Z,y + C,t + T), * = 1,2 (4.4) 

From this definition, one can generate the two-sided, three-dimensional, wavenumber-frequency 
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slope spectrum using the following relationship 


{ki, k 2 ,uj) = g-j J J J Rf exp [~iki£ — ik 2 ( — iur] d( dr, i = 1,2 

(4.5) 

where the wavenumber vector is k = (k cos 9, k sin 9) = {ki,k 2 ), and to is the wave frequency. 
The total slope spectrum is the sum of the two components, 

$ sl (k h k 2 ,u)) = $f(k u k 2 ,Lu) + $ s 2 l (k u k 2 ,u}) (4.6) 

The two-sided wavenumber spectrum can be obtained by simply integrating the three- 
dimensional spectrum in the frequency domain, 

fOO 

$ (M) = / $ sl (k,9,u>)duj (4.7) 

J — OO 

roo 

$ sl (k 1 ,k 2 ) = / $ sl (ki,k 2 ,Lj) dui (4.8) 

J—oo 

For further information on the physical design and data analysis development for the SLSG 
please reference [8] from which the above description was taken. 

The Fourier transform described in equation 4.5 cannot be calculated using straightfor¬ 
ward fast Fourier transform (FFT) techniques for the same reason that the similar equation 
for estimating the circular wire wave gage two-dimensional surface height wavenumber spec¬ 
trum cannot be: the data points in autocorrelation space are unevenly distributed. The 
circular wire wave gage processing overcomes this limitation by fitting a two-dimensional 
surface to the autocorrelation function and then estimating the spectrum. The processing 
algorithm used for the SLSG uses a different method. Since the SLSG data is sampled in 
evenly spaced increments in time, the temporal portion of the spectral integral is performed 
first using FFT methods. The resulting two-dimensional spatial integral is then numerically 
integrated using finite area elements associated with each lag coordinate. The finite areas 
are different for each spatial lag magnitude. A detailed description of this technique is given 
in the dissertation of Mete Uz [66] from the University of Rhode Island. 

Two sets of data were obtained from Dr. Hara with the help of one of his graduate 
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students, Hua Wei. The first data set is the fully processed three-dimensional slope spectrum 
in 6 minute segments from the CoOP ‘97 experiment. The second data set was composed 
of the calibrated time series slope data for the same 6 minute segments used to generate the 
first data set. To compare the processing algorithms, which will be referred to as SLSGproc 
and CWWGproc, the second data set was modified so as to allow it to be processed using 
the CWWGproc. The modifications included subsampling it spatially by using only every 
third point around the 15 centimeter diameter sampling circle used by the SLSG. The raw 
SLSG data is collected at 150 points around this circle, so the second data set was composed 
of only 50 sample points for each scan. This subsampling was necessary for two reasons, 
the first being to reduce the CWWGproc processing time which increases significantly as 
the number of points around the sampling circle becomes larger, and the second was to 
attempt to make the data set simulate something very similar to what the CWWG would 
have measured had it been actually sampling the identical wave field as the SLSG. The 
second data set was also bandlimited in frequency. This was done, again, to simulate 
CWWG conditions because during the CWWG calibration the measured wave heights are 
bandlimited before they are used to estimate the spectrum. 

SLSG data from Julian day 198 in 1997 was processed with the SLSGproc algorithms and 
the resulting three-dimensional spectrum $ sl (k x ,k y ,uj) was reduced to a two-dimensional 
spectrum by integrating out the frequency variable as explained in equation 4.8. The 
resulting slope spectrum representing a total of 12 minutes of sampling is shown in Figure 4- 
26. The same slope data was then processed using the CWWGproc algorithms after being 
filtered as described in the previous section. The resulting slope spectrum representing 
the same 12 minutes of data is shown in Figure 4-27. The wind velocity during this 12 
minute segment of data collected by the SLSG was approximately 8.7 m/s at 18 meters 
above the sea surface. When the total mean square slope (mss) was determined from each 
slope spectrum estimate by calculating the area underneath each spectrum, there was a 
significant discrepancy between the CWWGproc and SLSGproc results. The CWWGproc 
results gave a mss of 0.033 while the SLSGproc results gave a mss of 4.08 * 10 -6 . For the 
given wind speed during this data set of 8.3 m/s at 10 meters 3 the expected mean square 

3 Adjusted from the 18 meter wind speed using the Charnock model. 
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Figure 4 26: Scanning Laser Slope Gage Normalized Two-Dimensional Slope Spectrum 
Processed with SLSGproc from Data Sets ds023255.198 and ds023909.198. The spectrum 
is plotted as log 10 (4 >si ). The units are in (rad) 2 /(rad/cm) 2 . The magnitude of the 
spectrum was normalized as discussed in section 4.2. 
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Figure 4-27: Scanning Laser Slope Gage Two-Dimensional Slope Spectrum Processed with 
CWWGproc from Data Sets ds023255.198 and ds023909.198. The spectrum is plotted as 
log 10 ($ si ). The units are in (rad) 2 /(rad/cm) 2 . 
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slope should be between approximately 0.03 and 0.06 [30]. After a discussion with Dr. Bock 
it was determined that the SLSGproc results contained a calibration error which accounted 
for the much lower than expected mean square slope value. This calibration error affects 
only the magnitude of the spectrum and not the angular distribution. Since the primary 
purpose of this processing method comparison was to assess the ability of the CWWGproc 
to accurately reproduce the directional characteristics of a spectrum, the magnitude of 
the spectrum generated by the SLSGproc algorithms was normalized. This normalization 
was implemented by integrating out the angular variable in both the CWWGproc and 
SLSGproc slope spectra resulting in a one-dimensional (1-D) slope spectrum as a function 
of wavenumber for each method. The SLSGproc 1-D spectrum was then multiplied by a 
constant to make the area underneath it equal to the area underneath the CWWGproc 1-D 
spectrum. 

The two spectra shown in the logarithmic plots in Figures 4-26 and 4-27 are very sim¬ 
ilar (the logarithmic scale significantly enhances small differences). A closer look at their 
shapes is shown in Figure 4-28 which compares constant wavenumber slices from each of 
the processing methods. As can be seen, the CWWGproc at the lower wavenumbers does 
not appear to reproduce the narrowness of the spectrum which is evident in the SLSGproc 
results. It also seems to smooth the undulations visible in SLSGproc spectrum. This is ex¬ 
pected given that the input data for the CWWGproc algorithms is a subset of the SLSGproc 
input data (only every third data point is used) so the angular resolution of the CWWGproc 
results should be less than for the SLSGproc results. Additionally, the smoothing window 
used in the CWWGproc data has a width of about 0.4 rad/cm which would also account 
for some smoothing of the spectrum. 

The key purpose for comparing the shapes of the spectra using the two processing 
methods was to show that given the same input data, the two algorithms produced ap¬ 
proximately the same result. The above figures show this to be generally true. It will be 
assumed, therefore, that based on this comparison, the CWWGproc algorithm for generat¬ 
ing the two-dimensional spectrum is corroborated although the CWWGproc methods show 
a tendency to broaden the spectrum at low wavenumbers. 
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Figure 4-28: Directional Comparison Between SLSGproc and CWWGproc. The solid curves 
represent slices of the spectrum from the CWWGproc and the dashed curves represent 
slices of the spectrum from the SLSGproc. The two numbers next to each curve are the 
wavenumbers (in rad/cm) at which the slice is taken from the CWWGproc and SLSGproc 
spectra, respectively. The slices are not taken at the exact same wavenumbers because the 
spectral estimates are generated at a limited number of wavenumbers which are different 
for each processing method. 
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4.3 Discussion 


It is important to view the data presented in this chapter as a body of knowledge about 
wind and wave conditions at the Buzzards Bay Tower on October 28, 1998 (except for 
the last section about the scanning laser slope gage). Consistency in the data is therefore 
important. Evidence of this exists between the sonic anemometer measurements and the 
tower anemometer data: both show similar evolution of the wind direction over the course 
of the experiment from about 210° to 225°N. The scatterometer measurements also indicate 
a change in wind direction evidenced by the fact that the peak energy level shifts with the 
mean wind direction very consistently. Additionally, the long wave data measured by the 
capacitive wire wave gage indicates that the long wave energy is propagating primarily 
from the Southwest toward the Northeast. (It should be noted that in the meteorologic 
community wind direction is reported as the direction from which the wind is blowing, 
while in the surface wave community wave energy is most often reported as the direction 
in which the waves are propagating.) Finally, the single wire frequency spectra calculated 
from both the capacitive wire wave gage wires and the circular wire wave gage wires shows a 
remarkable match in both magnitude and slope over the ranges in which they both measure 
data as shown in Figure 4-29. This supports the accurate calibration of both instruments, 
as the frequency spectrum is a very reliable and consistent measure of ocean surface energy 
density. This also supports the assertion that the general statistical measures collected 
with the circular wire wave gage are the same as those collected at the tower even though 
they were not exactly co-located. The differences between the scatterometer directional 
measurements and the CWWG directional measurements will be discussed in the following 
chapter. 

How this body of data both with its consistencies and questions relates to the current 
state of knowledge in the fields of short surface wave theory and ocean backscatter modeling 
will be discussed in the following chapter. 
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Figure 4-29: Comparison of the Frequency Spectrum as Measured by the Capacitive Wire 
Wave Gage and the Circular Wire Wave Gage. The solid line is the frequency spectrum 
from the capacitive wire wave gage data set 10281414 (see Figure 4-13). The dots are the 
frequency spectrum from the circular wire wave gage data set 10281515 (see Figure 4-20). 
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Chapter 5 


The Wave Height Spectrum and 
Scatterometry 


the data presented in the last chapter offer a unique opportunity to study and compare 
in situ measurements of high-resolution two-dimensional wavenumber spectra and electro¬ 
magnetic backscatter return collected co-temporally. No currently published experiment or 
study has accomplished this. The results from this experiment are mixed. The scatterom- 
eter data compares well with current theories and other scatterometer measurements, but 
the circular wire wave gage two-dimensional spectrum results compare favorably with only 
some of the current spectral models. And with respect to the angular alignment of the wave 
gage energy distribution with the wind direction, the data is inconsistent with the expec¬ 
tation that the wind and short waves are generally well aligned. This chapter will focus on 
comparing the separate pieces of the Buzzards Bay Tower data with current models, and 
•will then treat the data as a collective unit to compare with models which attempt to use 
portions of the data (such as a two-dimensional wavenumber spectrum) as input to more 
complex models (such as electromagnetic scattering models); the goal being to explain as 
much of the deviation from the current theories as possible. 
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5.1 The Wave Height Spectrum 


Frequency and wavenumber representations of the surface wave height spectrum provide 
similar information about the energy distribution of surface waves, but they are collected 
and often modeled in different ways, especially in the gravity-capillary region of the spec¬ 
trum where Doppler shifting makes transformations between the two representations am¬ 
biguous. This section is therefore broken into two subsections, the frequency spectrum and 
the wavenumber spectrum, in order to show the evolution of current models in each area. 

5.1.1 Frequency Spectrum 

The omnidirectional frequency spectrum is one of the easiest statistical descriptions of 
surface wave height to determine because it can be calculated from a single time series of 
surface wave height measurements. Phillips [48] was the first to present a widely accepted 
model for the analytical form of this spectrum. He theorized that if the duration and fetch 
of the wind are large, then the smaller-scale components of the wind-generated wave field 
may be in a condition of equilibrium determined by the requirements for attachment of 
the crests of the waves. His model is based on a dimensional analysis of this idea of an 
equilibrium range and is valid for the rear face (high frequency portion) of the spectrum. 
It takes the form 


S(u) ~ ag 2 u 5 , u> peak < u; < — (5.1) 

where g is the acceleration due to gravity, u is frequency in rad/s, u* is the wind friction 
velocity in m/s, and a is a parameter independent of u and environmental conditions, which 
is obtained by fitting experimental data. 

While data at the time appeared to substantiate this theory, further study indicated 
that a actually depended in some manner on other parameters such as fetch and wind 
speed or friction velocity. Pierson and Moskowitz [49] extended the Phillip’s spectrum to 
cover the long wave portion of the spectrum and modeled it after their own data for fully 
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developed seas: 


S(w) = ag 2 uj 5 exp(-/3(u> 0 /u;) 4 ), tu 0 = g/U (5.2) 

where U is the wind speed at 19.5 meters above the sea surface and a and j3 are constants to 
be determined. The next major advance came with the analysis of data from the Joint North 
Sea Wave Project (JONSWAP) conducted in 1968 and 1969 which included fetch-limited 
seas. The resulting analysis prompted the modification of the Pierson-Moskowitz spectrum 
by multiplying it by a ‘peak enhancement’ factor q 6 [27] where 7 is a fitted parameter 
describing the degree of peakedness and 8 is an exponential dependent on frequency (a 
variable), peak frequency of the spectrum (a constant), and a fitting parameter 07 which 
describes the width of the peak region. 

It had been accepted by this time that the spectrum was a function of some combination 
of energy balancing factors such as input energy from wind, energy transfer between waves 
of different wavelength, and viscous dissipation. Most models after this time include some 
form of fetch or wave age dependence. Continuing analysis and at-sea experimentation in 
the 60’s and 70’s lead to a more fundamental modification of the basic Phillips spectrum, 
however. Much of this newer data supported a slope of lo~ 4 [62, 46, 6 ]. Donelan, Hamilton, 
and Hui [17] (DHH) provided the next wide-reaching and in-depth study of surface waves 
based on both lake (Lake Ontario) and laboratory (CCIW) experimentation in the late 70’s. 
Their conclusions have been widely accepted as the basis upon which current models are 
tested. Their study gave strong evidence to support the newer u>~ 4 power law and was used 
to generate a model which modified the JONSWAP spectrum to account for the different 
power law. It takes the form 

S{u) = ag 2 u)~ 4 u)~ 4 ak exp 

with 

T = exp [—(a; - uj peak ) 2 /2a 2 u 2 eak ] 

where the four parameters of the model are c o vea ^ the peak frequency of the spectrum, a, 


) 4 ' 


u 


7 


(5.3) 
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a fitted rear face variable, 7 , a fitted peak enhancement factor, and <7, a fitted peak width 
variable. Donelan et al. [17] also give approximate relationships between the fitted variables 
and the phase speed at the peak of the spectrum as well as the wind speed. 

An example of the frequency spectrum measured by the circular wire wave gage (CWWG) 
and the capacitive wire wave gage is compared with the Donelan, Hamilton, and Hui spec¬ 
trum in Figure 5-1. It is evident from this comparison that the measured spectrum near the 
spectral peak is closer to the DHH uj~ 4 , but at frequencies higher than about 1 Hz the slope 
increases to near tu~ 5 . It should be noted, however, that the DHH work only modeled the 
spectrum up to frequencies approximately three times the spectral peak and the extension 
out to much higher frequencies was done for illustrative purposes only. It stands to reason 



Figure 5-1: Comparison of the DHH Frequency Spectrum Model with a Frequency Spectrum 
Measured at the Buzzards Bay Tower. The solid line is the frequency spectrum from the 
capacitive wire wave gage data set 10281315 (see Figure 4-13). The dots are the frequency 
spectrum from the circular wire wave gage data set 10281415 (see Figure 4-24). The dashed 
line is the DHH spectral model given in equation 5.3. The dash-dot line is the DHH spectrum 
modified to use an w~ 5 rear face slope (similar to the original JONSWAP spectral model). 
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that the resulting difference between the measured spectrum and the DHH model simply 
indicates that at the higher frequencies, other factors besides wind speed and the peak phase 
speed (which are used as indicators of wave age) may be affecting the spectrum, or that the 
energy balances may be different. From a dimensional standpoint, if the slope of the mea¬ 
sured omnidirectional wavenumber spectra of -2.5 (shown in Figure 4-24) is transformed 
into frequency space using the dispersion relation and no Doppler shifting, the resulting 
frequency spectrum would have a slope of cj“ 4 , which indicates that some of the steepening 
in the frequency spectrum may be the result of Doppler shifting, although there are other 
physical mechanisms at these frequencies that could affect the spectral slope including wind 
drift and ambient current [17]. There have been convincing arguments and measurements 
in the literature, however, that support a frequency spectral slope of somewhere between 
-4 and -5, but the limited amount of wavenumber and meteorological data available with 
this data set do not allow for suggested modifications to or specific support of any current 
models except in the matter of the overall slope of the frequency spectrum. 

Most of the frequency spectrum models developed after the mid-80’s were not based on 
direct measurements of the frequency spectrum, but were produced using indirect methods 
involving transformations or integrations of either a one- or two-dimensional wavenumber 
spectrum or a more complex three-dimensional spectrum (4>(fc, cj)). These models are based 
on measurements or analytical development in wavenumber space, so they will be presented 
in the following section. 

5.1.2 One- and Two-Dimensional Wavenumber Spectra 

It has been more difficult to obtain experimental measurements of the wavenumber spectrum 
due to the need for extensive spatial simultaneous sampling of wave elevations. As a result, 
the earliest models of the wavenumber spectrum were generated by transforming a frequency 
spectrum model using some form of the dispersion relation. The problem with this method 
arises at high frequencies due to Doppler shifting of the short waves components by the long 
wave orbital velocities causing transformations to become ambiguous. Instrumentation has 
been developed over the last 10 or so years which has finally allowed the wavenumber 
spectrum at higher frequencies to be measured directly, which has lead to more robust 
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wavenumber spectrum models. The circular wire wave gage discussed in this thesis is one 
such instrument. 

There are two commonly used forms for describing the two-dimensional wavenumber 
spectrum, the direct form and the spreading function form. The direct form itself has two 
representations: $(k x ,ky) = <£>(fc) and 3>(fc, 9). The total mean square wave height energy, 
rj 2 , can be obtained from either form using the following equation 

/ oo roo rir ro o 

/ $(fc)dfc=/ / ${k,9)kdkd9 (5.4) 

-oo J —OO J—IT JO 

where k = (k x , k y ) is the wavenumber vector, k = |As] is the magnitude of the wavenumber 
vector, and 9 is the angle of the vector k. The omnidirectional spectrum shall be defined as 

S{k)de f = [* $(k,9)kd$ (5.5) 

describing the wave height spectral density as a function of wavenumber, k, for all directions. 
The spreading function form of representation is: 

9(k,$) = $(k,e max )D(k,9) ( 5 . 6 ) 

where D(k,9 ) is a non-dimensional spreading function used to model the azimuthal dis¬ 
tribution of energy and is conventionally presented with D(k, 0) = 1 for all k, 1 where the 
wind is assumed to be sufficiently steady so that the dominant wind waves are aligned with 
the mean wind direction (9 = 0). The function $(fc, 9 max ) will be called the ‘wavenumber 
spectral slice in the peak energy direction’ in this text and represents the distribution of 
energy purely as a function of wavenumber. A number of directional spreading functions, 
D(k,9 ) were described in detail in Chapter 2, but for thoroughness Tables 5.1 and 5.2 list 
these spreading functions again and include their associated high frequency wavenumber 
spectral slice models, $(k, 9 rnax ). 

The two-dimensional wave height spectra generated using the CWWG data repre¬ 
sent the directional distribution of energy in a purely spatial domain (lacking any fre- 

1 This definition is different from spreading functions which require a unity polar integral. 
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Year Reference 


Spreading Function 


1990 Banner [3] 


1994 Apel [ 1 ] 


D(k,9) = cosh ~ 2 [p{9 - 9 max )\ 

j 2.28{k/k p )-°- 65 > 0.97 < k/k p < 2.56 

10-0-4+0.839 exp[-0.567 In .(*/*,)], fc/fcp > 256 

D(k , 0) = exp{—(# - 0 up ^ nd ) 2 /[20 2 (/c)]} 
l/[20f(fc)j = 0.14 + 5.0(/c/fc p )- 1 - 3 


1996 Caudal et al. [9] D(k,9) = exp(-c*0 2 ) 

a = 0.14 + 5.0(fc/A:p)- 1 3 + 5, 0.9 < k < 1180 

£ _ lo(-0 029C/+0.064)X' 2 + (0.1271/+0.050)2f+(—0.126C/—0.930) 

X = lo §io( fc ) 

1997 Elfouhaily et al. [19]“ D(k , 0) = (27r)- 1 [l + A (k) cos(20)] 

A{k) = tanh[0.173 + 4(c/c p ) 2 - 5 + 0.565u*(0.23/c) 25 ] 

1997 Romeiser et al. [55] D(k, 6) = exp[-0 2 /(2<5 2 )] 

l/(2<5 2 ) = 0.14 + 0.5[1 - exp(—kU /400)]+ 

5exp[2.5 - 2.6 In (17) - 1.31n(ifc)] 


Table 5.1: Directional Spreading Functions from Various Models, k is the wavenumber in 
rad/m, k p is the wavenumber at the peak of the wave height spectrum in rad/m, U is the 
wind speed at 10 meters above the water surface in m/s, c is the phase speed in m/s, Cp is 
the phase speed at the peak of the wave height spectrum in m/s, 9 is the azimuthal direction 
in radians, and u* is the friction velocity in m/s. 

“The Elfouhaily spreading function model does not conform to the typical definition for the spreading 
function given in equation 5.6 in the sense that D(k, 0) ^ 1 for all k. Elfouhaily does choose, however, to 
compare his (non-normalized) spreading function and the resulting omnidirectional spectrum (as defined in 
equation 5.5) to other spreading functions and omnidirectional spectra in the literature which do require 
D(k, 0) = 0, so for the purposes of comparison herein, the model he presents in his paper will be used 
unchanged. 
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Year Reference Wavenumber Spectral Slice in the Peak Energy Direction 


1990 Banner [3] 4>(fc, 6 max ) = 0.45 x 10 4 (U/c p ) 1 / 2 k 4 , k in cycles/meter 

1994 Apel [1] 9 max ) — 0.00195 exp[—(k p /k) 2 }J p Hik~ 4 

j _ ^ yexp[-(v/fc--y/^p) 2 /(°- 32 fcp)] 

Hi = {[1 + (fc/100) 2 ]- 1 + exp{ln[—49.5 + 34.5(1 - exp(—U/4.7))]} 
0.8fccosh[(A; — 400)/450] _1 }exp(—fc 2 /6283 2 ) 


1997 Elfouhaily $(Mmas) = [5/ + B h ]fc“ 4 

et al. [19] B h = 0.00115(1 + 3 ln(ii70.23)]c _1 exp[-0.25(fc/371 - l) 2 ] 

Bi = 0.00 3^U/^p c p /cexp[-1.25(/c p /fc) 2 ] 

[1.7 + 61og(C//cp)] exp ^Iv / ^p- 1 ] 2 [° 16 + a64 (Wcp) -3 ]" 1 } 
exp[—0.316U/c p (^Jk/kp — 1)] 


1997 Romeiser $(Jfc, 0 max ) = P L {k, U)W H {k)U^ k h- 4 

et al. [55] P L = 0.000195exp{-(fc 2 /fc 2 ) + 0.53exp[-(vX- ^) 2 /{0.32k p )}} 

W H = [ 1 + (fc/280) 7 ' 2 ] 0,5 [l + (fc/75) 22 ] _1 [l + (fc/1300) 32 ] -2 
exp(— k 2 /8885 2 ) 

(3 = [1 — exp(—fc 2 /183 2 )] exp(—/c/183)+ 

[1 — exp(—A:/33)] exp{-[(fc — 140)/220] 2 } 


Table 5.2: High Frequency Wavenumber Spectral Slice in the Peak Energy Direction from 
Various Models, k is the wavenumber in rad/m except where noted, k p is the wavenumber 
at the peak of the wave height spectrum in rad/m, U is the wind speed at 10 meters above 
the water surface in m/s, c is the phase speed in m/s, Cp is the phase speed at the peak 
of the wave height spectrum in m/s, 6 is the azimuthal direction in radians, and u* is the 
friction velocity in m/s. 
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quency/temporal information); so the spectral energies in directions differing by 180° are 
the same. Guissard [22] considers this and explains that there are actually two possible 
interpretations of the notion of a spectrum which lead to different functions that do not 
have the same angular behavior: one is symmetric (the wave height spectrum), the other is 
not. He continues on to explain that for the purposes of electromagnetic scattering mod¬ 
els, although there are often upwind-downwind asymmetries in backscattered return, this 
difference cannot be attributed to the contribution of the spectrum through linear Bragg 
scattering mechanisms. This is because these mechanisms are directly determined by the 
wave height spectrum which is, as defined, symmetric. So the alternate (non-symmetric) 
form of the spectrum, while useful, is not appropriate for use in linear scattering models. 2 
For the purposes of consistency in the spectral comparisons between the CWWG data 
and the models, therefore, only the (symmetric) wave height spectrum will be used. The 
spreading models presented in Table 5.1 are not constrained to be symmetric, though, so a 
symmetrized or folded spreading function is used instead. This folded spreading function 
is defined as 


D* s (k, e ) = 0.5 [D(k, 6) + D(k, e - tt )3 (5.7) 

where the spreading functions D(k,6 ) are those given in Table 5.1. Finally, to permit 
effective comparison between the models and the CWWG data, since the resulting folded 
spectra of the previously non-symmetric models will have somewhat smaller magnitudes 
at D*(k, 0), 3 all of the spreading models and the CWWG azimuthal slice data have been 
normalized as a function of k by their upwind magnitude for consistency with the previous 
convention that D(k, 0) = 1. This normalized spreading function will be defined as 


D s (k,6) = D*{k,0)/max[D* s (k,e max )} (5.8) 


2 Guissard does offer some possible non-Bragg mechanisms to explain the upwind-downwind asymmetry 
in his paper. It should be noted that some very recent and more complicated scattering models than 
are presented in this thesis do attempt to incorporate a method to account for this asymmetry while still 
including the symmetric wave height spectrum dependencies required for linear Bragg scattering theory. 

3 This is true because the ‘upwind’ magnitude of D(k, 0) is constrained to be equal to 1, but the ’downwind’ 
magnitude of D(k, 0) is much less than 1 (i.e. the non-symmetry) so the sum of the ‘upwind’ and ‘downwind’ 
components multiplied by 1/2 will be less than 1. 
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where the maximum is calculated for each value of k and occurs when 8 rnax = 0 for each of 
the models and where 9 max = the maximum of the measured wave height spectrum for the 
circular wire wave gage data. 

Figures 5-2, 5-3, and 5-4 compare the omnidirectional spectrum and wavenumber spec¬ 
tral slice models with the data collected with the circular wire wave gage. The Banner 
model is not shown in either the two-dimensional spectral slice plot or the omnidirectional 
plot because it is nearly identical to the Apel model. The omnidirectional spectra calculated 
from the models were evaluated numerically using 

S(k) = r $(fc, 9 max )D s (k , 9) k d9 (5.9) 

J —7r 

A Caudal model spectrum in the dominant wind direction is not listed in the table or shown 
in the plots because Caudal et al. did not propose their own form of the spectrum, instead 
they suggested using Apel’s $(k, 9 max ), which is what was used to calculate the S(k ) integral 
which is shown for the Caudal model. Finally, the CWWG data has a resolution down to 
approximately 0.2 rad/cm due to instrument and processing constraints. The windowing 
used during the processing, therefore, causes the slight dip in the spectra near the low 
wavenumber cutoff. 

It is apparent from Figures 5-2 through 5-4 that the Romeiser et al. model comes 
closest to accurately describing the spectral values calculated from the circular wire wave 
gage data both in magnitude and slope. Interestingly, at wavenumbers above about 1 or 2 
rad/cm the CWWG spectral slice data begins to reduce its slope. This could be explained 
a number of ways. Banner [3] presented some analytical results indicating that the slope of 
the spectrum near the spectral peak can reduce the slope of the high frequency end of the 
spectrum through long wave spatial (wavelength) modulation of the short waves. He also 
showed that ambient surface currents can have the same effect. Thus, given the likelihood 
of strong currents in the location of the experiment, there is a possibility that the surface 
current velocities were larger than the phase velocities of the short waves, leading to what 
may be a significant effect on the expected spectral slope and energy distribution. Another 
explanation could be increased instrumentation error at these smaller amplitude waves. 
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Figure 5-2: Comparison of the Omnidirectional Spectrum Models and Wavenumber Spec¬ 
trum Slice Models with Data Collected at the Buzzards Bay Tower (10281415). The thick 
solid line is a slice of the two-dimensional wavenumber frequency spectrum in the max¬ 
imum spectral direction as measured with the circular wire wave gage (CWWG). The 
CWWG measurements are averaged over 10 minutes for the omnidirectional spectrum, 

and 20 minutes for the spectral slice. The thin lines represent the model functions: - 

= Apel,-= Elfouhaily et ah, and • • • = Romeiser et al. The thin solid line in the 

omnidirectional spectrum plot is the Caudal et al. model. The solid line with the small 
square data points is a reference curve for each plot showing a particular wavenumber slope. 
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Figure 5-3: Comparison of the Omnidirectional Spectrum Models and Wavenumber Spec¬ 
trum Slice Models with Data Collected at the Buzzards Bay Tower (10281445). The thick 
solid and thick dashed lines are slices of the two-dimensional wavenumber frequency spec¬ 
trum as measured with the circular wire wave gage (CWWG). The solid line (B) is a slice 
in the mean wind direction, and the dashed line (A) is a slice in the maximum spectral 
direction. The CWWG measurements are averaged over 10 minutes for the omnidirectional 
spectrum, and 20 minutes for the spectral slice. The thin lines represent the model func¬ 
tions: - = Apel,-= Elfouhaily et al., and • • • = Romeiser et al. The thin solid line 

in the omnidirectional spectrum plot is the Caudal et al. model. The solid line with the 
small square data points is a reference curve for each plot showing a particular wavenumber 
slope. 
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Figure 5-4: Comparison of the Omnidirectional Spectrum Models and Wavenumber Spec¬ 
trum Slice Models with Data Collected at the Buzzards Bay Tower (10281515). The thick 
solid and thick dashed lines are slices of the two-dimensional wavenumber frequency spec¬ 
trum as measured with the circular wire wave gage (CWWG). The solid line (B) is a slice 
in the mean wind direction, and the dashed line (A) is a slice in the maximum spectral 
direction. The CWWG measurements are averaged over 10 minutes for the omnidirectional 
spectrum, and 20 minutes for the spectral slice. The thin lines represent the model func¬ 
tions: - = Apel,-= Elfouhaily et al., and • • • = Romeiser et al. The thin solid line 

in the omnidirectional spectrum plot is the Caudal et al. model. The solid line with the 
small square data points is a reference curve for each plot showing a particular wavenumber 
slope. 
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The high wavenumber tail of the spectral slice in the maximum spectral direction (as well 
as those in the mean wind direction) shift between 20 minute data sets in a nonuniform 
manner, however, indicating that instrumentation noise is probably not entirely responsible 
for the variations. The Elfouhaily et al. model is the only other model shown which 
comes close to reproducing the spectral slopes of the at sea data. Interestingly, and not yet 
interpreted, is that the Romeiser model is based on fitting a model to radar backscatter 
data, while the Elfouhaily model is based on wave tank measurements. The Apel model 
is based on both scatterometer and tank data, but is valid only for fully developed seas, 
so it is not unexpected that it generates too steep a spectral slope. Likewise, the Caudal 
model is also too steep, and while it was based on scatterometer data, it is not limited to 
only fully developed seas. Another reason for the slope inconsistencies between the CWWG 
data and the models in the spectral slice graphs may be the sensitivity that the models have 
to the location of the spectral peak. The frequency of the spectral peak at the Buzzards 
Bay Tower, found using the capacitive wire wave gage, may differ from the true frequency 
of the spectral peak by up to 5% because the location of the spectral peak is a random 
variable which is generally under sampled (especially during limited duration experiments 
such as this one) as was explained by Donelan and Pierson in 1983 [16]. The key information 
in these plots is that the CWWG omnidirectional spectra support the —2.5 to —3 slopes 
suggested in the current models, but that the CWWG wavenumber spectral slices support 
a slope of —3 which is smaller than the —4 model slopes. 

Figures 4-21, 4-22, and 4-23 from Chapter 4 show that the wind and the peak energy 
of the measured wave height spectra are not always well aligned. To see this in detail, 
Figures 5-5, 5-6, and 5-7 show slices of the two-dimensional wavenumber spectrum at a 
constant wavenumber (resulting in a circular slice of the 2-D wavenumber spectrum) for both 
the non-dimensional normalized form of the spreading models, D s (k, 6), and the normalized 
CWWG data. In all of these plots the peaks of the spreading models are centered on the 
mean wind direction, which is their standard orientation, while the data is presented as 
measured. 

These plots show a significant shifting of the peak location over the course of the experi¬ 
ment. Although the central measure studied in this thesis is the spectral width, a digression 
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Figure 5-5: Comparison of the Normalized Directional Spreading Models with Data Col¬ 
lected at the Buzzards Bay Tower (10281415). The thick solid line in each plot is a constant 
wavenumber slice (a circular ring) of the two-dimensional wavenumber frequency spectrum 
measured by the circular wire wave gage (CWWG). The CWWG measurements are aver¬ 
aged over 20 minutes. The wavenumber above each plot is the wavenumber at which the 

slice was taken. The thin lines represent the model functions:-= Elfouhaily et al., • • • = 

Romeiser et al., and — = Caudal et al. The Apel spreading model is not included because 
it is effectively isotropic at these wavenumbers when symmetrized using equation 5.7. The 
plots are angularly centered near the maximum experimental data peak, but the model 
directional slices are centered on the mean wind direction which is shown as a vertical 
dashed line at 213°. All the models and the CWWG data are normalized as described in 
equation 5.8 
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Figure 5-6: Comparison of the Normalized Directional Spreading Models with Data Col¬ 
lected at the Buzzards Bay Tower (10281445). The thick solid line in each plot is a constant 
wavenumber slice (a circular ring) of the two-dimensional wavenumber frequency spectrum 
measured by the circular wire wave gage (CWWG). The CWWG measurements are aver¬ 
aged over 20 minutes. The wavenumber above each plot is the wavenumber at which the 

slice was taken. The thin lines represent the model functions:-= Elfouhaily et al., • • • = 

Romeiser et al., and — = Caudal et al. The Apel spreading model is not included because 
it is effectively isotropic at these wavenumbers when symmetrized using equation 5.7. The 
plots are angularly centered near the maximum experimental data peak, but the model 
directional slices are centered on the mean wind direction which is shown as a vertical 
dashed line at 221°. All the models and the CWWG data are normalized as described in 
equation 5.8 
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Figure 5-7: Comparison of the Normalized Directional Spreading Models with Data Col¬ 
lected at the Buzzards Bay Tower (10281515). The thick solid line in each plot is a constant 
wavenumber slice (a circular ring) of the two-dimensional wavenumber frequency spectrum 
measured by the circular wire wave gage (CWWG). The CWWG measurements are aver¬ 
aged over 20 minutes. The wavenumber above each plot is the wavenumber at which the 

slice was taken. The thin lines represent the model functions:-= Elfouhaily et al., • ■ • = 

Romeiser et ah, and — = Caudal et ah The Apel spreading model is not included because 
it is effectively isotropic at these wavenumbers when symmetrized using equation 5.7. The 
plots are angularly centered near the maximum experimental data peak, but the model 
directional slices are centered on the mean wind direction which is shown as a vertical 
dashed line at 226°. All the models and the CWWG data are normalized as described in 
equation 5.8 
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is presented here to consider possible explanations for the shift in energy direction. There 
are a number of possible explanations for this energy shift: one, the compass data was not 
incorporated properly, and two, the meteorological conditions at the location of the buoy 
shifted over time and were not represented well by measurements made at the Buzzards 
Bay Tower. With respect to the first, the compass data collected on the buoy was used to 
reference the local measurements, which shifted as the buoy floated in the waves, to a geo¬ 
graphic coordinate system. It was further verified that this referencing was done properly, 
and the actual orientation of the buoy over the course of the experiment was assessed and is 
shown in Figure 5-8. These plots show that while the buoy did rotate over the course of the 
experiment, it did not progressively rotate further and further in one direction which could 
have affected the continual shift of the peak energy location through possible errors in the 
processing. With respect to the second possible explanation, it is curious that the first data 
set is fairly well aligned with the wind measurements made at the tower, but that as the 
buoy moved further from the tower, the following two data sets show a systematic rotation 
of the peak energy. The winds further from the tower may actually have shifted, but it is 
difficult to believe that within three nautical miles of the tower they would have shifted by 
up to 90°. Another possible explanation was mentioned earlier as a possible explanation for 
the unexpectedly low wavenumber slope in the wavenumber slice plots - surface currents. 
These currents could have affected the directional properties of the short waves, but it is 
difficult to explain the long term shifting unless the buoy was caught in some kind of quasi¬ 
stable eddy of sorts. Unfortunately, there is no way to know for sure as no surface wind or 
surface current measurements were collected at or near the buoy. Given that the CWWG 
instrument and processing errors which could have caused these spurious results have been 
rigorously checked, the precise reason for this behavior is currently unknown. 

Figures 5-5, 5-6, and 5-7 illustrate the variability in the peak location of the CWWG 
data well, but it is difficult to determine if the widths of the measured spectra and the 
models are comparable, so a calculation of spectral width as a function of wavenumber was 
generated. A mean square angular width was defined as 


Width(k ) = 


j(e-e CG ) 2 Hk,e)de 

f$(k, 6) dd 



(5.10) 
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(a) 



(b) 



(c) 


Figure 5-8: Buoy Orientation, (a) is from data set 10281415, (b) is from data set 10281445, 
and (c) is from data set 10281515. The angles given are the direction that the compass was 
facing, but the compass was not aligned with the wind vane on the buoy, so these angles 
do not correspond to the wind directions presented elsewhere in this thesis. 
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Figure 5-9: Comparison of the Widths of the Directional Spreading Models with Data 
Collected at the Buzzards Bay Tower (10281415). * = CWWG data, + = Banner model, 
o = Romeiser model, < = Caudal model, o = Elfouhaily model The widths are shown in 
degrees. 20 minutes of data were used to generate these width measurements. 


where 

(5.11) 

and T*(/c, 6) is the two-dimensional wavenumber spectrum. 

Figure 5-9 shows a calculation of the widths of the CWWG data and the models for 
data set 10281415. This Figure clearly illustrates that using this measure of width, the 
CWWG data is narrower than the models predict except at the lowest wavenumber bin. 
Specifically, the Romeiser and Caudal models propose a slight narrowing of the spectrum 
as a function of wavenumber, while the Banner and Elfouhaily models do not. The circular 
wire wave gage data, thus, supports an even greater narrowing of the spectrum than the 
Romeiser and Caudal models. 

Figures 5-5, 5-6, and 5-7 show some abberations in the measured spectra at the higher 
wavenumbers (they are not very smooth), though. To verify that the width measurements 


je$(k,e) de 

f$(k,0) d9 
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have converged and would not simply get wider if more data were averaged, a measure of 
width as a function of averaging time was generated and is shown in Figure 5-10. An issue 
that is raised by the plots in Figure 5-10 is the significant (up to 50°) variability of the width 
measurements for short averaging times. Inspection of the two-dimensional wavenumber 
spectra for these short time period averages indicates that the large widths are not generally 
indicative of actually broad spectra, but rather are primarily generated by multi peaked 
spectra which have energy distributed in two distinct regions. It is possible that a different 
measure of width, say a 3-dB estimator may show a different kind of variability, but some 
method would still need to be developed to deal with the multi peaked distribution. 

5.2 Electromagnetic Scattering Models 

The scatterometer data collected at the tower offers a slightly different opportunity to see 
how well the collected data agrees with models which use the surface wave height spectra 
as an input parameter. The data collected by the X-band scatterometer at the tower 
exhibits similar azimuthal distribution as other scatterometer data sets available. Figure 5- 
11 shows the Buzzards Bay Tower scatterometer azimuthal measurements on a decibel 
scale [dB = 10 * log 10 (power)] for the 10281315 data set. The upwind-to-crosswind ratio 
is approximately 6 for the HH polarization and between 5 and 6 for the VV polarization. 
Similar ratios exist for the other two scatterometer data sets at the tower. These upwind- 
to-crosswind ratios are nearly identical to those found by Masuko et al. [44] for similar 
meteorological conditions and incidence angles. Additionally, Feindt et al. [20] include a 
plot (their Figure 16) of upwind-to-downwind ratios as a function of wavenumber which 
accurately predicts the ratios measured by the scatterometer at the tower, however they 
show a large scatter in the data. Interestingly, even accounting for the scatter in the data, 
the Feindt plot shows an increasing upwind-to-crosswind a 0 ratio as a function of Bragg 
wavenumber, indicating a narrowing trend in the scatterometer azimuthal data at higher 
scatterometer frequencies. 

A detailed description of the form of the electromagnetic scattering model relevant to the 
experimental measurements presented herein was given in Chapter 2. A rather generic form 
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(c) 10 minute averages 


(d) 20 minute averages 


Figure 5-10: Convergence of Circular Wire Wave Gage Width Measurements. The width 
values are calculated using all three spectral data sets (10281415, 10281445, and 10281515). 
















Figure 5-11: Directional X-Band Scatterometer Return from Data Set 10281315. This plot 
is a replica of Figure 4-17 except that it is shown on a decibel scale. The upper curve and 
data points are VV polarization and the lower fitted curves and data points are for HH 
polarization. The power has units of I 2 +Q 2 where I and Q are measured in volts. The solid 
line represents a best fit of the data to the function cos 2s while the dashed line 

represents a best fit of the data to the function cosh 2 — @max))- The s and (3 values for 
the HH and VV polarizations are given in the Figure 4-17 caption. The data points shown 
are 85 point averages. 
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of this model for linear wave theory which includes the now well established Bragg scattering 
foundation along with a tilt modulation factor is repeated here with the normalized radar 
backscatter cross section (NRCS) being given as 

0o(<M) = 871-fco J G(0,7Lx,7L2/)r(^ ,Q)P{7Lx,lLy)dlLxd'Uy (5.12) 

where (p is the incidence angle, 6 is the azimuthal look direction, ko is the electromagnetic 
wavenumber, G{<j>,iLx->lLy) is a non-dimensional polarization dependent geometric reflec¬ 
tion coefficient, T(<p', 6) is the folded surface displacement wavenumber spectrum evaluated 
at the Bragg wavelength, and -P( 7 Lx> 7 Ly) is the probability density of the long wave slopes, 
7 L x and 7 l v - The only two azimuthal dependencies given in this equation are within the 
folded wavenumber spectrum and the long wave probability distribution. 

Using Lyzenga’s [41] simplification of Plant’s [50] Gaussian approximation of the slope 
probability distribution detailed in Chapter 2, it was concluded that a skewed probability 
distribution does have the ability to impact the angular dependence of the NRCS. Using 
extrapolated data from the CWWG two-dimensional wavenumber spectrum for the shape 
and magnitude of the short wave input to the scattering model, and using the long wave 
measurements from the capacitive wire wave gage to estimate the long wave slopes, the 
output of the modeled scatterometer return with the tilting is narrower than the simple 
Bragg scattering model return without tilting effects, but the effect is minimal. 

Based on the simplified scattering model in equation 5.12 then, the primary cause of the 
breadth of the scattered return is the width of the wave height spectrum. This is supported 
by the Buzzards Bay Tower data as shown in Figures 5-12 and 5-13. Both figures show 
width calculations for both the scatterometer data, the corresponding CWWG data, and 
the previously described models. As can be seen, the width of the scatterometer data is 
on the same scale as the widths measured by the CWWG data, but it is still narrower 
than current spectral models predict. It can be concluded, therefore, that the circular wire 
wave gage data shows support for the current scattering theories based primarily on Bragg 
scattering. 
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Figure 5-12: Comparison of the Widths Between Spreading Models, the Circular Wire 
Wave Gage, and the X-Band Scatterometer (10281415(CWWG) and 10281315(scat)). * = 
CWWG data, + = Banner model, o = Romeiser model, « = Caudal model, o = Elfouhaily 
model, t> = HH scatterometer return width, and x = W scatterometer return width. The 
scatterometer widths are plotted at their approximate Bragg wavenumber. The widths are 
shown in degrees. 
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Wavenumber (rad/cm) 


Figure 5-13: Comparison of the Widths Between Spreading Models, the Circular Wire 
Wave Gage, and the X-Band Scatterometer (10281445(CWWG) and 10281344(scat))). * = 
CWWG data, + = Banner model, o = Romeiser model, < = Caudal model, o = Elfouhaily 
model, t> = HH scatterometer return width, and x = VV scatterometer return width. The 
scatterometer widths are plotted at their approximate Bragg wavenumber. The widths are 
shown in degrees. 
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5.3 Conclusions 


The one- and two-dimensional wavenumber spectral data presented in this and the previous 
chapter agree with only some of the current model predictions. Specifically, the measured 
one-dimensional frequency spectra show a slope which increases from to near w -5 at the 
higher wavenumbers. This agrees with the DHH spectral model at the low frequency end, 
but behaves much more like earlier models predicting the steeper slope of —5. However, 
experimental data in the literature at these frequencies is not consistent and it is believed 
that the slope of the frequency spectrum lies somewhere between -4 and -5 depending 
on environmental conditions. Meanwhile, the omnidirectional wavenumber spectrum plots 
show the data supporting the —2.5 to —3 slopes which are typical in three of the models 
shown. But in the wavenumber spectral slice plots the CWWG data shows a gradually 
decreasing slope which is probably associated with physical processes, such as surface cur¬ 
rents, which are not incorporated into the models presented here. The wavenumber spectral 
slice models, thus, predict a slope of near -4, but the CWWG data never gets below about 
-3. The cause for the differences between the measured and the modeled frequency and 
wavenumber slopes are not well understood, but it is believed that both the orbital veloc¬ 
ities of the long waves and possible strong surface currents may be somewhat responsible. 
There is not enough evidence within these data sets to argue this conclusively, however, 
and further investigation of the surface currents at the Buzzards Bay Tower may shed some 
light on these discrepancies. 

The measured width of the spectra show a narrowing as a function of wavenumber which 
is supported, to a slightly lesser degree, by only two of the models presented herein. The 
other models discussed do not support any narrowing as a function of wavenumber. Addi¬ 
tionally, the X-band scatterometer data measured at the tower compare very well to other 
available experimental scatterometer data sets. It was also shown that a simplified compos¬ 
ite surface model based on Bragg scattering could reproduce the width of the normalized 
radar cross-section as measured by the X-band scatterometer at the Buzzards Bay Tower, 
when given as input an interpolated form of the two-dimensional wave height spectrum. 

The only other in situ data avaliable at these wavenumbers was published by Hwang et 
al. in 1996 [31] and Hara et al. in 1998 [24]. Both Hwang and Hara used scanning laser 
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slope gages. These instruments measured slope instead of wave height, and both were able 
to calculate what is known as a saturation spectrum. The saturation spectrum is gener¬ 
ally defined a s k A * S(k) where S(k) is the omnidirectional wave height spectrum, or as 
k 2 *S sl (k) where S' sl (k) is the omnidirectional wave slope spectrum. This spectrum is often 
used instead of the wave height spectrum for comparative purposes because if it is assumed 
that there is a saturation limit in the slope of the energy distribution of k 4 as described in 
Chapter 2, then the saturation spectrum should have zero slope. The saturation spectra 
from these two papers at wind friction velocities similar to those experienced during the 
Buzzards Bay Tower experiment are about two times smaller than the CWWG saturation 
spectra, but they show similar wavenumber slopes. Field data collected by these authors 
show large scatter in magnitude of the spectrum, however, and it is not believed that the 
Buzzards Bay Tower omnidirectional spectrum data constitute any significant divergence 
from the previous experimental data, especially since the other surface conditions for these 
experiments are not well matched to those experienced at the Buzzards Bay Tower. Given 
that these data sets were collected in differing ocean conditions it isn’t appropriate, un¬ 
fortunately, to compare the directional distribution of energy from them with the CWWG 
data. 

Incidentally, a concern about the interpretation of the data with respect to the current 
models exists because of the shifting of the peak energy direction away from the mean wind 
direction. As mentioned previously, this may be due to some currently unknown error in the 
instrument, a true shift in the mean wind direction near the buoy, or a significant surface 
current effect, but there is nothing to either support or refute any one of these explanations, 
given that there were no local wind or surface current measurements made at the buoy. The 
shift may also indicate that at these short wavelengths there is the possibility of significant 
amounts of energy propagating in directions other than the mean wind direction, but this 
is hard to explain over the course of a 20 minute average. It is important to note that 
circular wire wave gage measures the short wave distribution over a patch of the sea surface 
which is small with respect to the long waves. The resulting averaged spectra represent an 
estimate of the variable energy distribution over a very small spatial domain. The statistical 
variability of this process at these scales is not well understood, and without additional and 
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more lengthy studies it is difficult to draw any firm conclusions about why the averaged 
energy distribution is not well aligned with wind measurements in the synoptic area. 

While the data presented in this thesis is not comprehensive enough to make conclusive 
recommendations for slight modification to some of the current models to account for the 
small variations which are shown in both magnitude of the spectrum and its width, it 
does indicate that further studies need to be conducted. Specifically, further experimental 
data with both radars and in situ instruments capable of measuring the two-dimensional 
wavenumber spectrum over a significant amount of time should be performed. Only through 
more thorough investigation will the discrepancies between the models and data of this 
nature finally be resolved. 
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Chapter 6 


Summary 


Short waves on the ocean surface play an important role is surface electromagnetic (e.m.) 
scattering. Electromagnetic scattering theory is a well developed field, but the short wave¬ 
length portion of the surface wave field has only recently been experimentally explored, and 
a single, consistent model of the wave height spectrum has yet to be developed. A new 
instrument was developed to measure the height of waves with 2-30 cm wavelengths at an 
array of locations which can be post-processed to generate an estimate of the surface wave 
height spectrum. This instrument was deployed in an experiment which gathered not only 
in situ measurements of the two-dimensional wave height spectrum, but also coincident 
scatterometer measurements. This thesis describes the development of this instrument and 
the comparison of the data which was collected during this experiment with current spectral 
and electromagnetic surface theory models. 

This work began with a description of surface wave theory including a basic description 
of surface wave propagation and current theories specifying the two-dimensional energy 
distribution out to the gravity-capillary region of the spectrum. A presentation of electro¬ 
magnetic scattering theory was then given. A brief historical outline of its development was 
given as well as the mathematical basis to understand the current models. This chapter 
ended with an account of other instruments available which measure surface wave spectra. 

Chapter 3 introduced the circular wire wave gage. It talked about the buoy structure 
and housings, the layout, modification, and purpose of the analog electronics, the data 
acquisition hardware, the basic operation of the instrument, as well as the post-deployment 
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data processing necessary to generate the wave height spectrum estimate. 

The last two chapters described the Buzzards Bay Tower experiment and the comparison 
if the resulting data with current models. The instrumentation usage and layout description 
was followed with plots and images of all of the meteorological, buoy, and scatterometer 
data collected. The circular wire wave gage data processing was also compared with that 
of another directional spectrum measuring instrument to validate its methods. The wave 
height spectrum as measured by the circular wire wave gage was then compared to the 
models presented in Chapter 2. It was found that while the omnidirectional spectrum 
measurements were similar in slope but slightly larger in magnitude than the models, the 
spectral slice from the in situ data in the dominant wave direction showed variability both in 
slope and magnitude. This was evident also in the two-dimensional wavenumber spectrum 
estimates generated from the CWWG data. It was shown that the azimuthal spreading 
of these spectral plots was only slightly underestimated by the current spectral models. 
The data collected at the Buzzards Bay Tower was also useful for comparison with current 
scatterometer models. Specifically, it was shown that the directional spreading measured 
by the circular wire wave gage when extrapolated out to the Bragg wavelength of the X- 
band scatterometer used, was very near the corresponding azimuthal width measured by 
the scatterometer. Given that the azimuthal dependence in most current linear theory 
composite surface theory scattering models is entirely derived from the spreading of the 
surface wave height spectrum, this data presents new evidence that these models are very 
close to accurately modeling the Bragg spreading mechanisms. 

6.1 Future Work 

The circular wire wave gage has provided new information about short wave distribution on 
a platform that is easy to deploy. There are certainly improvements which would enhance 
its performance and usefulness for short wave studies, however. These include adding 
a ‘motion-pack’ system capable of monitoring pitch, heave, and roll, allowing long wave 
information to be ascertained which could potentially provide the ability to subsample the 
data at just the wave peaks or troughs allowing actual in situ measurements of the differences 
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between the short surface wave energy distribution over the long wave profile. A wind sensor 
mounted on the buoy would also help with temporal correlation of the wave data with the 
wind data. Additionally, surface currents and long wave orbital velocities may affect the 
short wave spectra as was suggested in Chapter 5, so future deployments incorporating 
colocated surface current measurements would be very informative. Finally, while two- 
dimensional wavenumber spectra provide very useful data about energy distribution on 
the ocean surface, this information could be significantly enhanced with the addition of 
information in the third dimension: time. It is possible to obtain the three-dimensional 
wave height spectrum from the circular wire wave gage data up to a frequency of nearly 50 
Hz given that each scan is collected at 100 Hz, but analysis and data processing algorithms 
have not yet been developed for this purpose. 

While the Buzzards Bay Tower experiment generated new information about the validity 
of current wavenumber spectrum and electromagnetic scattering models, it did not provide 
enough information to generate model corrections or to suggest exactly which processes 
caused the differences between the experimental data and the models. Additional studies 
are needed to look at the directional characteristics of short sea waves over a series of varying 
meteorological conditions in conjunction with azimuthal scattering, specifically because the 
variability in directional properties of short sea waves may be significant and is highly 
dependent on wind and wave conditions. 
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Appendix A 


Spectral Model References 

Table A.l: Chronological Listing of Spectral Model References (1950’s - 1990’s), 
exp. = experimental, spread. = azimuthal spreading, spec. = spectrum 
Most of the articles references here are from one of the following journals: Journal of 
Geophysical Research, Journal of Physical Oceanography, Journal of Marine Research, or 
the Journal of Fluid Mechanics. Full reference information for some of these articles can be 
found in the bibliography. 


Author (s) 

Development 

Output Form 

Phillips [1958] 

Cox [1958] 

dimensional analysis 

u spec 

Pierson & Moskowitz [1964] 
Ewing [1969] 

Valenzuela et al. [1971] 

exp. (Moskowitz data) 

uj spec. 

Black [1974] 

theoretical 

co spec. 

Tyler et al. [1974] 

Stillwell et al. [1974] 

exp. (SAR fit) 

2D k spec. 

Mitsuyasu et al. [1975] 

exp. (cloverleaf buoy) 

2D a; & spread. 

Kitaigorodskii et al. [1975] 

theoretical 

uj & k spec. 


continued on next page 
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continued from previous page 


Author(s) 

Development 

Output Form 

Larson & Wright [1975] 

Long &; Huang [1976] 

Plant & Wright [1977] 

exp. (tank/e.m. backset.) 


Wu [1977] 

Regier et al. [1977] 

Mitsuyasu [1977] 

Reece, Jr. [1978] 


2D slope spec 

Bjerkaas & Riedel [1979] 

Lleonart et al. [1980] 

theoretical & experimental 

u & k spec. 

Hasselman et al. [1980] 

Kahama [1981] 

exp. (JONSWAP/pitch-roll buoy) 

2D w & spread. 

Huang et al. [1981] 

exp. (JONSWAP) 

k spec. 

Fung & Lee [1982] 

experimental (Radar) 

u> & k spec. 

Kitaigorodskii [1983] 

theoretical & experimental 


Holthuijsen [1983] 

exp. (N. Sea/stereophotography) 

2D k & spread. 

Tang & Shemdin [1983] 

exp. (MARSEN/laser slope follow) 

u spec. 

Komen et al. [1984] 

Jackson et al. [1985b] 

theoretical 

2D u spec. 

Hwang [1985] 

exp. (Atlantic/laser slope) 

2D k spectrum 

Phillips [1985] 

theoretical & experimental 

2D k spectrum 

Donelan et al. [1985] 

exp. (Lake Ontario/wave staff fit) 

uj spec., spread. 

Haimback & Wu [1986] 

exp. (Lab/scanning slope gage) 

prob. dist. spread. 

Donelan & Pierson [1987] 

exp. 


Atakturk et al. [1987] 

exp. (Lake WA./wave gage) 

to spec 

Shemdin & Hwang [1988] 

exp. (MARSEN & TOWARD/laser) 

slope, uj spec. 

Shemdin et al. [1988] 

exp. (TOWARD/stereophoto.) 

2D k sp'ec. 


continued on next page 
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continued from previous page 


Author(s) 

Development 

Output Form 

Banner et al. [1989] 

exp. (BHP platform/stereophoto) 

2D k spec 

Jahne [1989] 

exp. (Delft tank/opt grad k si) 

slope u k k spec. 

Banner [1990] 

exp. (previous data) 

ID spec, spread. 

Jahne k Riemer [1990] 

exp. (Delft tank/optical grad.) 

2D k slope spec. 

Bock k Hara [1992] 



Klinke k Jahne [1992] 



Jackson et al. [1992] 

exp. (radar) 

mean sq. slope 

Apel [1994] 

exp. (Donelan, Banner, Jahne) 

k spec. 

Caudal k Proud’hom [1994] 

theoretical dvlpmnt of <j 0 

k spec, anisotropy 

Young [1994] 

theoretical dvlpmnt of array config. 

spread. 

Hara et al. [1994] 

exp. comparison btw gage & radar 

tu-k slope spec. 

Hara & Plant [1994] 

exp. (SAXON k MARSEN/radar) 

MTF 

Zhang [1995] 

exp. (Scripps tank/optical grad.) 

k spec 

Young et al. [1995] 

exp. (Lk. George, AU wave gage) 

direc. spread. 

Liu k Yan [1995] 

theoretical dvlpmnt/ERS-1 comp. 

k spec, k spread. 

Caudal k Hauser [1996] 

exp. (JONSWAP, Donelan) 

direc. spread. 

Kitaigorodskii [1996] 



Hwang et al. [1996] 

exp. (HiRes ARI) 

k spec. 

Elfouhaily et al. [1997] 

exp. (JONSWAP, Phillips, Kitai.) 

2-D k spec. 

Hara et al. [1997] 

exp. (Canada flume/laser slope) 

w-k si spec., spread. 

Hwang [1997] 

exp. (Atlantic/laser slope, radar) 

k spec. 

Romeiser et al. [1997] 

exp. (airborne scat.) 

k spec., spread. 

Hanson et al. [1997] 

exp. (catamaran/laser slope) 

2D uj spec. 
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Appendix B 


Scattering Model References 

Table B.l: Chronological Listing of Scattering Model References (1950’s - 1990’s), 
exp. = experimental, spread. = azimuthal spreading, spec. = spectrum, theo. = theoretical 
Most of the articles references here are from one of the following journals: Journal of Geo¬ 
physical Research, Journal of Physical Oceanography, Journal of Marine Research, IEEE 
Transactions on Geoscience and Remote Sensing, or the Journal of Fluid Mechanics. Full 
reference information for some of these articles can be found in the bibliography. 


Author(s) 

Development 

Frequency(s) 

Rice [1951] 

Saxton et al. [1952] 

theoretical (perturbation method) 


Wright [1966] 

experimental 

X-band 

Barrick et al. [1968] 

theoretical (composite surf, theory) 


Wright [1968] 

theoretical & exp. (tank) 

P,L,C,X-band 

Bass et al. [1968] 

Valenzuela & Laing [1970] 

Wright & Keller [1971] 

Long [1974] 

Klein & Swift [1977] 

theoretical (perturbation method) 



continued on next page 
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Author (s) 

Development 

Frequency(s) 

Brown [1978] 



Valenzuela [1978] 

Jones et al. [1978] 

Alpers & Hasselmann [1978] 
Alpers et al. [1981] 

Fung & Lee [1982] 

full theo. development 

400-8910 MHz 

Plant & Keller [1983] 

Lyzenga et al. [1983] 

Wentz et al. [1984] 

Kwoh & Lake [1984] 

Keller et al. [1985] 

Hasselmann et al. [1985] 

exp. (MARSEN), MTF 

X,L-band 

Durden & Vesecky [1985] 

Plant [1986] 

Masuko et al. [1986] 

Holliday et al. [1986] 

Fung & Pan [1987] 

Pan & Fung [1987] 

exp. (RADSCAT &; others) 


Donelan & Pierson Jr. [1987] 
Phillips [1988] 

Wetzel [1990] 

Brown [1990] 

full theo. development/data sigma-o 

L,C,X,Ka band 

Plant [1990] 

Trizna et al. [1991] 

exp. (ocean, lake/Doppler spec.) 

1.5-14.0 GHz 

Giovanangeli et al. [1991] 
Jessup et al. [1991a] 

exp. (azimuthal tank studies) 

36 GHz 


continued on next page 
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Author(s) 

Development 

Frequency(s) 

Jessup et al. [1991b] 

Kasilingam et al. [1992] 

Ebuchi et al. [1993] 

Wetzel [1993] 

Trizna et al. [1993] 

Apel [1994] 

exp. (lab) 

5.5-35 GHz 

Hara & Plant [1994] 

Liu & Pierson [1994] 

Plant et al. [1994a] 

exp. (SAXON & MARSEN), MTF 


Lee et al. [1995] 

exp. (Scotland/Doppler spec.) 

X-band 

Trizna [1996] 


X-band 

Trizna [1997] 

theoretical (Brewster angle) 

L,S,C,X,Ku,Ka-band 

Romeiser & Alpers [1997] 

theoretical 


Plant [1997] 

exp. (SAXON), Doppler 

X-band 
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Appendix C 


Wave Measuring Instrument 
References 


Table C.l: Chronological Listing of Wave Measuring Instrument References (1950’s-1990’s). 
exp. = experimental, spread. = azimuthal spreading, spec. = spectrum 
The Output Form column also indicates whether the work is based primarily on field ex¬ 
periments or whether it is based on laboratory work. Most of the articles references here 
are from one of the following journals: Journal of Geophysical Research, Journal of Phys¬ 
ical Oceanography, Journal of Marine Research, or the Journal of Fluid Mechanics. Full 
reference information for some of these articles can be found in the bibliography. 


Author(s) 

Instrument Type 

Output Form 

Cox & Munk [1954] 

Aerial Photogrammetry 

field, slope stats. 

Cote et al. [1960] 

Aerial Photogrammetry 


Barber [1961] 

Array Directional Resolution 

theoretical development 

Longuet-Higgins et al. [1963] 

Floating Buoy 

field 

Prettyman et al. [1969] 

Laser Slope Gage 


Stillwell [1969] 

Wu et al. [1969] 

Single Image Photography 



continued on next page 
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continued from previous page 


Author(s) 

Instrument Type 

Output Form 

Garret [1969] 

Capacitive Probe 

field, f spec. 

McGoldrick [1970] 

Capacitive Probe 


Barrick [1972] 

High Frequency (HF) Radar 


Tober et al. [1973] 

Laser Slope Gage 

lab, l o slope spec. 

Sturm & Sorrell [1973] 

Laser Slope Gage / 1 pt. 

lab 

Sugimori [1975] 

Stereo Photography 


Hughes et al. [1977] 

Laser Slope Gage 


Palm et al. [1977] 

Laser Slope Gage 


Davis k Regier, L. [1977] 

Wire Array Design (resistive) 

field, cu-k spec. 

Schuler, D. [1978] 

Dual Frequency Radar 

field 

Lubard et al. [1980] 

Opt. grad. vs. Laser Slope 


Wu et al. [1981] 

Laser Scanning Slope Gage 


Lange et al. [1982] 

Ampl. vs. Slope Comparison 

lab, to slope spec. 

Tang k Shemdin [1983] 

Optical Slope Gage (on wave follower) 

field, f spec. 

Keller k Gotwols [1983] 

Optical 


Holthuijsen [1983] 

Stereo Photography 


Donelan et al. [1985] 

Capacitive Wire Wave Array 

lab/field, 2D k spec. 

Haimbach k Wu [1985] 

Optical Scanner 


Hasselman et al. [1985] 

Synthetic Aperture Radar 


Walsh et al. [1985] 

Surface Contour Radar 


Ebuchi et al. [1987] 

Specular Scattering 

lab 

Shemdin et al. [1988] 

Stereophotography 

field, 2D k spec. 

Smith k McCardle [1990] 

Capacitive Multiwire Buoy 

field, 2D k spec. 

Jahne k Riemer [1990] 

Surface Light Gradient 

lab 

Lyzenga et al. [1991] 

2-ant. SAR 


Li et al. [1993] 

Laser Slope Gage 

lab/field, 2D k spec 


continued on next page 
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continued from previous page 


Author(s) 

Instrument Type 

Output Form 

Hwang et al. [1993] 

Laser Slope Gage 

lab, 2D k spec 

Chapman k Monaldo [1993] 

Capacitive Wire Wave Gage 

lab/field 

Goldstein et al. [1994] 

SAR/Acoustic Doppler 

field, 2D oj spec. 

Zhang k Cox [1994] 

Surface Gradient Detector 

lab, 2D k si. spec. 

Sletten et al. [1994] 

Polarimetric Radar 

lab 

Bock k Hara [1995] 

Laser Slope Gage 

lab/field, w-k spec. 
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Appendix D 


Analog Electronics Interface Board 
Electrical Diagram 


The diagram shown in Figure D-l represents the electrical schematic for the circular wire 
wave gage analog electronics interface board. This board is located in the central analog 
electronics housing and acts as the interface between the wave wires, the analog electronics, 
and the data acquisition system located in a separate housing. 
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Figure D-l: Analog Electronics Interface Board 
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and Computation Lab 













































































































































Appendix E 


Oscillator Circuit Electrical 
Diagram 


The diagram shown in Figure E-l represents the electrical schematic for the circular wire 
wave gage oscillator. This small board is mounted on the larger wiring interface board in 
the central analog electronics housing, and its output is the buffered high-frequency carrier 
signal used by all the channels. This signal is referred to as the 10 kHz sine wave input on 
the other circular wire wave gage analog circuit diagrams. 
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Voltage Regulators 



Figure E-l: Oscillator Circuit Diagram 
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Appendix F 


Data Acquisition Design Criteria 


The data acquisition system was designed to meet the following criteria 

• Sample a set of data channels at 100 Hz. 

• Collect (scan) each set of 64 channels as quickly as possible (needs to be collected 
quickly enough that ’stationarity’ of the scan can be assumed) 

• Store collected data from a typical day’s deployment on board the buoy 

• Provide means to download the data from the buoy to another computer without 
having to disassemble the data acquisition system 

• Sample the on-board compass at least once every second 

• Sample the temperature probe at least once every 5 to 8 seconds 

• Sample the battery voltage regularly and shut off if below a specified cutoff level 

• Minimize power consumption so as to minimize the number of batteries needed 

• Acquisition system must fit inside an 8-inch diameter aluminium housing 

• Provide easy-to-use command line operation of data acquisition program which allows 
for calibration and for user to define data set length, interval between data sets, and 
absolute start time 
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Appendix G 


Data Acquisition Hardware 
Schematic 


The diagram shown in Figure G-l represents the general schematic for the circular wire wave 
gage data acquisition system. The PC104 form-factor motherboard and associated hardware 
are mounted to vertical struts which are secured to the endcap of the data acquisition 
aluminum housing. 
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E -O 



Figure G-l: Data Acquisition Hardware Schematic 
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Appendix H 


Data Acquisition Software 


The following two programs control the general data acquisition of the circular wire wave 
gage. The first program is wave . c which manages the actual analog to digital conversion 
and data storage. The second program, setup, c, directs the watchdog timer to turn the 
instrument on and off based on user input. 


H.l Wave.c 


/******************************************************************** 

* Copyright (c) 1997, Woods Hole Oceanographic Institution 

itc*********************************************************.^.^.^.^^^ 

/* 

* name: wave.c 

* date: 16 February 1997 

* auth: k. prada 

* 

* revision history: 

* 12 may 97 - revise power control for J2 use with relay 

* 

* This program performs data acquisition for the WAVE GAGE system. 

* It samples 64 channels of wave gage information using 4 A/D boards 

* It also samples a KVH compass unit. The A/D boards are all 

* DM5406-1 from Real Time Devices. The compass unit is sampled 

* serial using COM2. 

* 

* Analog data are sampled using the burst mode on the A/D board and DMA 
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* to fill blocks of main memory until a sequence is complete. 

* The data are double buffered and written to disk. This version 

* uses a default sequence interval of 15 minutes. 

* 

* Data statistics: 


* 


* 

Sample Rate: 

100 Hz 


* 

Channels/sample/board: 

16 

(64) 

* 

Bytes/sample/board: 

32 

(128) 

* 

Bytes/second/board: 

3200 

(12800) 

* 

Bytes/block/board: 

16384 

(131072) (5.12 seconds/block) 


* 

* command line options: 

* -b blocks number of blocks (32 = default) 

* -g gain amplifier gain (1 = default) 

* -i interval restart interval in minutes (15 = default) 

* -n start sampling now 

* -q quiet mode 

* -s disable shutdown after sampling 

* 

*/ 

#include <stdio.h> 

#include <stdlib.h> 

#include <stdarg.h> 

#include <string.h> 

#include < ctype.h> 

#include <dos.h> 

#include <conio.h> 

#include "wave.h" 

#include "com_m.h" 

#include "pwrctl.h" 

#include "rtd.h" 

struct header 

{ 

char title[33]; 
char date[9]; 
char time[9]; 
char batt [7] ; 
char fill [6]; 

> hd; 
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dtm st,t; 

unsigned char far *buff,*bf; 
unsigned char kvhl [120],kvh2 [120]; 
extern int Debug_Flag; 
extern int com_chars; 
static int Interval; 
static int stm,sth; 
static int Shutdown; 

static FILE *fo,*fkvh,*ftmp; 

#define BATT_FACT0R 4.0 

#define BATT.MINIMUM 10.0 

#define DEFAULT.INTERVAL 15 

#define DEFAULT.BLOCKS 32 


/* time structures */ 

/* data buffer pointer */ 

/* compass buffers */ 

/* verbose flag */ 

/* C0M1 compass character count */ 
/* restart interval in minutes */ 

/* shutdown flag */ 


/* battery binary to voltage factor */ 
/* minimum allowable operation volts */ 
/* every 15 minutes */ 

/* 16 DMA blocks = 10.24 minutes */ 


/*- ^i 

/* SYSTEM SHUTDOWN */ 

/*-*/ 


void stop(int code) 

{ 

stop_watchdog(); /* disable watchdog */ 

set_ds_watchdog(0); 
printf("\n\n"); 


exit(code); /* back to system */ 

> 

/* */ 

/* SET TO RESTART AT SPECIFIED TIME */ 

/*-*/ 


void restart(void) 
{ 

int i; 


if(sth >= 0) 
st.hr = sth; 
if(stm >= 0) 
st.min = stm; 

set_ds_alarm(&st); /* set the alarm for next wakeup */ 

debugO'Waiting 5 seconds for Keyboard Hit\n"); 
debugO'Press any key to Inhibit System Shutdown\n"); 
for(i = 0; i < 5; i++) 

{ 

if(kbhit()) 

{ 
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getchO ; 
stop(1); 

> 

debug("\r 0 /.2d ", (5 - i)); 
delay(1000); 


> 

debug ( 11 \n"); 

system_off(); /* put system to sleep */ 

> 

/*-*/ 

/* SYSTEM SHUTDOWN AND SCHEDULE NEXT SEQUENCE START TIME */ 

/*-*/ 


void quit(int code) 

{ 

int i,sets,cset; 


debug("\n\nAll Power Bits Off\n"); 

power(Oxff,OFF); /* disable all power */ 


if(!Shutdown) 
stop(code); 


/* if shutdown disabled */ 


sets = get_ds(40); 
cset = get_ds(41); 

++cset; 

put_ds(41,cset); 

debug ("°/,d of */.d Data Sets 

if(cset >= sets) 

{ 


/* get number of sets to do */ 
/* get number of sets done */ 
/* increment sets done */ 

/* restore sets done */ 
Complete\n",cset,sets); 

/* if all sets are done */ 


debug("Shutting Down W/0 Wakeup\n\n"); 
stop_alarm(); 
goto getout; 


debugO'Prior Start Time: %02d:'/,02d:%2d\n" ,st .hr,st .min,st. sec) ; 
st.min += Interval; /* add wakeup interval in minutes */ 

while(st.min >= 60) /* adjust to proper time */ 


st.hr++; 
st.min -= 60; 

> 

if(st.hr > 23) 
st.hr -= 24; 
debugC'Next Start Time: 
set_ds_alarm(&st); 


/* increment hour */ 

/* decrement minutes */ 

/* now adjust hour */ 

7,02d’/,02d\n" , st. hr, st. min) ; 

/* set the alarm for next wakeup */ 
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getout: 

debugC'Waiting 5 seconds for Keyboard Hit\n"); 
debugC'Press any key to Inhibit System Shutdown\n"); 
for(i =0; i < 5; i++) 

{ 

if(kbhit()) 

{ 

getchO ; 
stop(code); 

} 

debug("\r°/,2d ",(5 - i)); 
delay(lOOO); 

> 

debug(" \n"); 

system_off(); /* put system to sleep */ 


} 

/* - */ 

/* ABORT THE OPERATION */ 


void getout(void) 

{ 

rtdl_trigger(OFF); 
rtdl_dma_stop(); 
rtdl_reset(); 

rtd2_trigger(OFF); 
rtd2_dma_stop(); 
rtd2_reset(); 

rtd3_trigger(0FF); 
rtd3_dma_stop(); 
rtd3_reset(); 

rtd4_trigger(OFF); 
rtd4_dma_stop(); 
rtd4_reset(); 

fclose(fo); 
fclose(fkvh); 
fclose(ftmp); 
quit(1); 


} 

/* - */ 

/* WAIT FOR DMA DONE */ 

/* - */ 


/* disable A/D */ 

/* disable the DMA channel 
/* clear the board */ 


*/ 


/* disable A/D */ 

/* disable the DMA channel 
/* clear the board */ 


*/ 


/* disable A/D */ 

/* disable the DMA channel 
/* clear the board */ 


*/ 


/* disable A/D */ 

/* disable the DMA channel 
/* clear the board */ 


* / 
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void dma_wait(void) 

{ 

set_ds_watchdog(99); 

while(rtdl_dma_done() >= 0); 
rtdl_trigger(OFF); 
rtdl_dma_stop(); 
rtdl_reset(); 

while(rtd2_dma_done() >= 0); 
rtd2_trigger(0FF); 
rtd2_dma_stop(); 
rtd2_reset(); 

while(rtd3_dma_done() >= 0); 
rtd3_trigger(0FF); 
rtd3_dma_stop(); 
rtd3_reset(); 

while(rtd4_dma_done() >= 0); 
rtd4_trigger(OFF); 
rtd4_dma_stop(); 
rtd4_reset(); 


/* wait for DMA done for RTD1 */ 
/* disable A/D */ 

/* disable the DMA channel */ 

/* clear the board */ 

/* wait for DMA done for RTD2 */ 
/* disable A/D */ 

/* disable the DMA channel */ 

/* clear the board */ 

/* wait for DMA done for RTD3 */ 
/* disable A/D */ 

/* disable the DMA channel */ 

/* clear the board */ 

/* wait for DMA done for RTD4 */ 
/* disable A/D */ 

/* disable the DMA channel */ 

/* clear the board */ 


/* - 

/* WRITE ADC DATA TO DISK 

/* - 

void disk_write(int page,unsigned blk,unsigned bytes,int pack) 

{ 

unsigned segment,i,s,d; 


*/ 

*/ 

*/ 


set_ds_watchdog(99); 

segment = page * 0x1000; /* set segment address */ 

buff = MK_FP(segment,blk); 

if(pack) 

{ 

bf = MK_FP(segment,blk); 
s = 2; 
d = 1; 

ford = 0; i < bytes; i++) 


> 


bf[d] = bf [s] ; 
s += 2; 
d++; 
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> 

if(fwrite(buff,1,bytes,fo) != bytes) 

{ 

error ("\nFile Write Error For Block */,d @ %p\n",blk,buff); 
getoutO ; 

> 


> 

/* */ 

/* WRITE KVH COMPASS DATA TO DISK */ 

/* - */ 


void kvh_write(char *b,int n) 

if(fwrite(b,l,n,fkvh) != n) 

{ 

error("\nFile Write Error For KVH\n"); 
getoutO; 

} 


> 

/* - */ 

/* GET AND WRITE TEMPERATURE TO DISK */ 

/* - */ 


void temperature(void) 

{ 

float t; 

t = (float)ltc(O) * 0.00121; 
fprintf (ftmp, "7,6.3f\n" ,t); 


> 

/* - */ 

/* MAIN PROGRAM */ 

/* - +/ 


void main(int argc,char **argv) 

{ 

int c,i,j,gain,block,nblks,page,flag,segment,kvh_cnt; 
unsigned int battery; 

char fo_name[32],kvh_name[32],tmp_name[32],path[40]; 
float vbat; 

extern char *optarg; 
extern int optind; 


/* - */ 

/* ANNOUNCE THE PROGRAM */ 

- */ 


printf ("WAVE GAGE SYSTEM - Version '/.s 7.s\n" ,__DATE__ ,__TIME__) ; 
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/*-*/ 

/* ENABLE THE WATCHDOG */ 

/* - */ 

set_ds_watchdog(99); 
start.watchdog(); 

/* - */ 

/* PRESET SOME VARIABLES AND GET COMMAND LINE OPTIONS */ 

/* - */ 


get_ds_time(&st); /* get and save start date/time */ 

nblks = DEFAULT_BLOCKS; /* default number of DMA blocks */ 

gain =1; /* preset gain = 1 */ 

Interval = DEFAULT_INTERVAL; /* default interval */ 

Debug_Flag =1; /* default is verbose */ 

Shutdown =1; /* default is shutdown after mission */ 

stm = sth = -1; 

while((i = getopt(argc,argv,"g:i:b:h:m:qs?")) != -1) 

{ /* command line args */ 

switch(i) 

{ 

case ’g’: /* specify gain */ 

gain = atoi(optarg); 
debug ("Set Gain: °/ 0 d\n" ,gain) ; 
break; 

case ’i’: /* specify operational interval */ 

Interval = atoi(optarg); 

debugC'Set Restart Interval: 7,d minutes\n" ,Interval); 
break; 

case J b’: /* number of blocks */ 

nblks = atoi(optarg); 
debugC'Set Blocks: 7,d\n" ,nblks) ; 
break; 

case 'q’: /* specify quiet mode */ 

Debug_Flag = 0; 
debugC'Set Quiet Mode\n"); 
break; 

case ’s': /* disable shutdown */ 

Shutdown = 0; 

debugC'Disable Shutdown\n"); 
break; 

case ’h’: 

sth = atoi(optarg); 

debugC'Set Start Hour: ’/,d\n" ,sth) ; 

break; 

case ’m’: 
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stm = atoi(optarg); 

debugO'Set Start Minute: */,d\n" ,stm); 

break; 

case : /* help */ 

printf("-g set adc gain\n"); 
printf("-i set restart interval (minutes)\n"); 
printf("-b number of DMA blocks (5.12 sec/block)\n"); 
printf("-q set quiet mode\n"); 
printf("-h set start sample Hour\n"); 
printf("-m set start sample Minute\n"); 
printf("-s disable shutdown at end of interval\n"); 
printfO 1 -? this help printout\n\n"); 
return; 

> 

> 


/* - */ 

/* ALLOW USER TO BAIL OUT IF NECESSARY */ 

/* - */ 


printf("Press ESC within 5 seconds to AB0RT\n"); 
for(i =0; i < 5; i++) 

{ 

printf ("\r°/,2d ",(5-i)>; 
if(kbhit ()) 

{ 

if (getchO == 27) 

{ 

stop_watchdog(); 
set_ds_watchdog(0) ; 
stop(O); 

> 

> 

delay(1000); 

> 


/* - */ 

/* POWER ON - INITIALIZE THE A/D BOARDS AND SERIAL I/O */ 

/*-*/ 


debug("\rEnable Power\n"); 

power(J2,0N); /* analog section on */ 

delay(1000); /* short wait */ 

rtdl_init(); /* setup the A/D board */ 

rtdl.gain(gain); /* set gain */ 

rtdl_pacer_rate((long)25000); /* set burst pacer rate */ 

rtdl_trigger_rate((long)100); /* set 100 Hz trigger rate */ 
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rtd2_init(); 
rtd2_gain(gain); 
rtd2_pacer_rate((long)25000); 
rtd2_trigger_rate((long)100); 
rtd3_init(); 
rtd3_gain(gain); 

rtd3_pacer_rate((long)25000); 
rtd3_trigger_rate((long)100); 
rtd4_init(); 
rtd4_gain(gain); 
rtd4_pacer_rate((long)25000); 
rtd4_trigger_rate((long)100); 

kvh_init(); 


/* setup the A/D board */ 

/* set gain */ 

/* set burst pacer rate */ 

/* set 100 Hz trigger rate */ 
/* setup the A/D board */ 

/* set gain */ 

/* set burst pacer rate */ 

/* set 100 Hz trigger rate */ 

/* setup the A/D board */ 

/* set gain */ 

/* set burst pacer rate */ 

/* set 100 Hz trigger rate */ 

/* init the KVH compass port */ 


-- 

/* READ BATTERY VOLTAGE - ABORT IF BELOW MINIMUM SPEC 

/* - 

debug("Read Battery: "); 
battery = 0; 

for(i 0; i < 10; i++) /* do 10 sample average */ 

battery += ltc(l); /* g e t battery value */ 

battery /= 10; 

vbat = (float)battery * 00.00121 * BATT_FACT0R; 
debug(" °/,6.2f Volts\n" ,vbat) ; 

if(vbat < BATT.MINIMUM) /* if less than allowable */ 


-*/ 


*/ 


-*/ 


error("%6.2f Battery Voltage TOO LOW!\n",vbat); 
error("ABORT Mission With No Further Wakeup\n"); 
error(" In 5 Seconds\n\n"); 

for(i = 0; i < 5; i++) 
delay(1000); 

stop_alarm(); /* no wake-up */ 

if(Shutdown) 

system.off(); /* pu t system to sleep */ 

else 


stop(l) ; 

> 

-- 

/* CHECK TIME FOR NORMAL 

-- 

if(stm >= 0 || sth >= 0) 
restart(); 


START TIME OR NOW 


debug ("Sequence Start @ %02d:‘/ 0 02d: %02d\n" 

-- 


st.hr,st.min,st.sec) 


-*/ 


*/ 


-*/ 


-*/ 
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/* 

/*— 


OPEN FILES AND WRITE HEADER 


*/ 


- */ 


debug ("Data Capture: 7.d Blocks ”,nblks); 

sprintf (fo_name, ” , / t 02d , / l 02d7.02d , /.02d. WAV" , st.mon, st. day, st.hr, st. min) 
debug ("To File: 7,s\n" ,f o_name) ; 

fo = fopen(fo_name,"wb"); /* open the file for writing */ 

if(fo == NULL) 

{ 


error("Error Opening Output 
quit(l); 


File 7.s\n" ,fo_name) ; 

/* abort until next sequence */ 


sprintf (kvh_name, "7,02d7.02d7.02d7.02d. KVH", st. mon, st. day, st. hr, st. min); 
fkvh = fopen(kvh_name,"wb"); /* open the file for writing */ 

if(fkvh == NULL) 


/* 

/* 

/* 


errorO'Error Opening Output File %s\n",kvh_name); 
quit(l); /* abort until next sequence */ 


> 

sprintf (tmp_name, "*/ t 02d7.02d7.02d7.02d. TMP" , st. mon, st. day, st. hr, st. mxn) 
ftmp = fopen(tmp_name,"wb"); /* open the file for writing */ 

if(ftmp == NULL) 

{ 

errorO'Error Opening Output File 7«s\n" ,tmp_name) ; 

quit(1); /* abort until next sequence */ 

> 

get_ds_time(&t); /* start of sample time */ 

sprintf(hd.title,"WAVE GAGE DATA"); 

sprintf (hd. date, "7.02d/7.02d/7.02d" ,t.yr,t .mon,t .day) ; 

sprintf (hd. time,"7.02d: 7,02d: 7.02d" ,t. hr, t. min, t. sec) ; 

sprintf (hd.batt, "7.6.2f " ,vbat); 

sprintf(hd.fill," \r\n"); 

fwrite(&hd,l,64,fo); /* write the header */ 

sprintf(hd.title,"WAVE GAGE COMPASS"); 

fwrite(&hd,l,64,fkvh); /* write the header */ 

sprintf(hd.title,"WAVE GAGE TEMPERATURE"); 
fwrite(&hd,1,64,ftmp); /* write the header */ 


ACQUIRE DATA 


*/ 

*/ 

■*/ 


flag = 0; 
block = 0; 
while(l) 

{ 

rtdl_scan(16); rtdl_mux(0); 
rtd2_scan(16); rtd2_mux(0); 


/* preset buffer flag empty */ 

/* preset block number */ 

/* each block is 5.12 seconds */ 
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rtd3_scan(16); rtd3_mux(0); 
rtd4_scan(16); rtd4_mux(0); 
rtdl_dma_start(9,0,16384); 
rtd2_dma_start(7,0,16384); 
rtd3_dma_start(8,0,16383); 
rtd4_dma_start(6,0,16383); 

com_int_rcv(2,kvhl,(long)100,0); /* start compass capture */ 
if(flag) /* flag says a buffer is full */ 

{ /* write hi block buffers */ 

disk_write(9,0x8000,(unsigned)16384,0); 
disk_write(7,0x8000,(unsigned)16384,0); 
disk_write(8,0x8000,(unsigned)16384,1); 
disk_write(6,0x8000,(unsigned)16384,1); 
kvh_write(kvh2,(100-kvh_cnt)); 

} 

temperature(); 

dma_wait(); /* wait for DMA done */ 

kvh_cnt = com_int_cnt(2); /* get count remainder */ 

com_int_off(2); /* stop compass capture */ 

flag = 1; 

if(++block == nblks) /* if all blocks are done */ 

{ /* write final lo block */ 

disk_write(9,0,(unsigned)16384,0); 
disk_write(7,0,(unsigned)16384,0); 
disk_write(8,0,(unsigned)16384,1); 
disk_write(6,0,(unsigned)16384,1); 
kvh_write(kvhl,(100-kvh_cnt)); 
break; 

> 

rtdl_scan(16); rtdl_mux(0); 
rtd2_scan(16); rtd2_mux(0); 
rtd3_scan(16); rtd3_mux(0); 
rtd4_scan(16); rtd4_mux(0); 
rtdl_dma_start(9,1,16384); 
rtd2_dma_start(7,1,16384 ); 
rtd3_dma_start(8,1,16383); 
rtd4_dma_start(6,1,16383); 

com_int_rcv(2,kvh2,(long)100,0); /* start compass capture */ 

disk_write(9,0,(unsigned)16384,0); 

disk_write(7,0,(unsigned)16384,0); 

disk_write(8,0,(unsigned)16384,1); 

disk_write(6,0,(unsigned)16384,1); 

kvh_write(kvhl,(100-kvh_cnt)); 

temperature(); 

dma_wait(); /* wait for DMA done */ 
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/* get count remainder */ 

/* compass capture off */ 

/* if all blocks are done */ 


kvh.cnt = com_int_cnt(2); 
com_int_off(2); 
if(++block == nblks) 

{ 

disk_write(9,0x8000,(unsigned)16384,0) 
disk.write(7,0x8000,(unsigned)16384, 0) 
disk_write(8,0x8000,(unsigned)16384, 1) 
disk.write(6,0x8000,(unsigned)16384,1) 
kvh_write(kvh2,(100-kvh_cnt)); 
break; 

> 

} 

fclose(fo); 
fclose(fkvh); 
fclose(ftmp); 
quit (0); 


H.2 Startup.c 


/******************************************************************** 

* Copyright (c) 1997, Woods Hole Oceanographic Institution 

/* 

* name: setup.c 

* date: 12 may 1997 

* auth: k. prada 

* 

* revision history: 

* 

* Setup runtime parameters for the wave program 

* 

*/ 


#include 

#include 

#include 

#include 

#include 


<stdio,h> 
<stdlib,h> 
<time.h> 
<dos,h> 
<conio.h> 


#include <pwrctl.h> 
#include <wave.h> 
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dtm t; 

void main(void) 

{ 

int i,hour,minute,sets; 
clrscrO ; 

printf("\n\nWARNING!!!!\n\nThis program performs a system setup"); 
printf("\nThen it shuts the system down with alarm set!"); 
printf("\n\nPress <RETURN> to continue or <ESC> to quit...."); 
i = getchO ; 
if(i == 27) 

return; 

get_ds_time(&t); 

debug("\n\nCurrent Time: 7 t 02d:'/,02d:‘/,02d\n" ,t.hr,t.min,t.sec) ; 

printf ("Enter Start Hour.? "); 

scanf("%d",&hour); 

printf ("Enter Start Minute.? "); 

scanf ("°/,d" ,&minute); 

printf("Enter Number of Datasets..? "); 
scanf ("7 0 d" ,&sets) ; 

put_ds(40,sets); /* number of data sets to do */ 

put_ds(41,0); /* data set accumulator */ 

t.hr = hour; 
t.min = minute; 

set_ds_alarm(&t); /* set alarm */ 

printf ("Next Start Time: %02d7,02d\n" ,t .hr ,t .min) ; 
printf("Shutdown in 2 seconds "); 
for(i =0; i < 8; i++) 

{ 

printf("."); 
delay(250); 

> 

system_off(); 

> 
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Appendix I 


Capacitive Wire Calibration 


The capacitive wire wave gage was calibrated using simple linear calibration curves which 
were generated by connecting a known capacitance to each of the electronics modules and 
measuring their voltage outputs the results of which are shown in Table 1.1. Straight lines 
were fit to the calibration curves generating offsets and slopes for each wire. Additionally, 
when the full 63” of wave wire 1 was submerged in saltwater, the output voltage was 8.35 
volts. This conversion factor of length/voltage was used to generate a maximum capacitance 
value for each wire. It was assumed that the maximum capacitance for each wire was 
identical because the wires were manufactured at the same time and were all cut to the 
same length. The calibration is accomplished by multiplying the measured output voltage 
from each capacitive wire by the calibration slope, adding the offset, and then multiplying 
by the conversion factor of 19.17 cm/uF. 


Capacitance, 5% 

Wire 1 

Wire 2 

Wire 3 

Wire 4 

Wire 5 

Wire 6 

0 (uF) 

0.02 (volts) 

0.01 

0.02 

0.01 

0.01 

0.02 

1 

1.02 

1.01 

1.02 

1.01 

1.02 

1.01 

2 

2.02 

2.02 

2.03 

2.02 

2.04 

2.03 

3 

2.96 

2.95 

2.97 

2.96 

2.96 

2.97 

4 

4.03 

4.02 

4.04 

4.03 

4.02 

4.04 

5 

5.01 

5.00 

5.02 

5.01 

5.04 

5.02 


Table 1.1: Calibration Data for Capacative Wire Wave Gage 
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